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(54) Computer-aided visualization and analysis system for nucleic acid sequence evaluation 

(57) A computer system (1) for analyzing nucleic 
acid sequences is provided. The computer system is 
used to perform multiple methods for determining 
unknown bases by analyzing the fluorescence intensi- 
ties of hybridized nucleic acid probes. The results of indi- 
vidual experiments may be improved by processing 
nucleic acid sequences together. Comparative analysis 
of multiple experiments is also provided by displaying ref- 
erence sequences in one area (814) and sample 
sequences in another area (81 6) on a display device (3). 
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Description 
COPYRIGHT NOTICE 

A portion of the disclosure of this patent document contains material which is subject to copyright protection. The 
copyright owner has no objection to the xeroxographic reproduction by anyone of the patent document or the patent 
disclosure in exactly the form it appears in the Patent and Trademark Office patent file or records, but otherwise reserves 
ail copyright rights whatsoever. 

GOVERNMENT RIGHTS NOTICE 

Portions of the material In this specification arose in the course of or under contract nos. 92ER81 275 (SBIR) between 
Affymetrix. Inc. and the Department of Energy and/or H60081 3-1 , -2 between Affymetrix, Inc. and the National Institutes 
of Health. 

BACKGROUND OF THE INVENTION 

The present invention relates to the field of computer systems. More specifically, the present Invention relates to 
computer systems for visualizing biological sequences, as well as for evaluating and comparing biological sequences. 

Devices and computer systems for forming and using arrays of materials on a substrate are known. For example. 
PCT applications WO92/10588 and 95/11995, incorporated herein by reference for all purposes, describe techniques 
for sequencing or sequence checking nucleic acids and other materials. Arrays, fer pertorjnirig. these operations may be 
formed in arrays according to the methods of, for example, the pioneering techniques disclosed in U.& Patent Nos. 
5,445,934 and 5384,261, and U.8. Patent Application No. 08/249,188. each incorporated herein by reference for all 
purposes. 

According to one aspect of the techniques descrtoed therein, an array of nudeic acid probes is fabricated at known 
locations on a chip or substrate. A labeled nucleic acid is then brought into contact with the chip and a scanner generates 
an image file (also called a cell file) indicating the locations where the labeled nucleic acids bound to the chip. Based 
upon the image file and identities of the probes at specific locations, it becomes possible to extract information such as 
the monomer sequence of DNA or RNA. Such systems have been used to form, for example, arrays of DNA that may 
be used to study and detect mutations relevant to cystic fibrosis, the P53 gene (relevant to certain cancers). HIV. and 
other genetic characteristics. 

Improved computer systems and methods are needed to evaluate, analyze, and process the vast amount of Infor- 
mation now used and made available by these pioneering technologies. 

SUMMARY OF THE INVENTION 

An Improved computer-aided system for visualizing and determining the sequence of nucleic acids is disclosed. 
The computer system provides, among other things, improved methods of analyzing fluorescent image files of a chip 
containing hybridized nucleic acid probes In order to call bases In sample nucleic add sequences. 

According to one aspect of the invention, a computer system is used to identify an unknown base in a sample nucleic 
acid sequence by the steps of: 

- inputting multiple probe intensities, each of the probe intensities being associated with a nucleic acid probe; 

- the computer system comparing the multiple probe Intensities where each of the probe intensities Is substantially 
proportional to a nucleic acid probe hybridizing with at least one nucleic acid sequence; and 

calling the unknown base according to the results of the comparison of the multiple probe Intensities. 

According to one specific aspect of the invention, a higher probe intensity is compared to a lower probe intensity to 
caD the unknown basa According to another specific aspect of the invention, probe Intensities of a sample sequence 
are corrpared to probe intensities of a reference sequence. According to yet another specific aspect of the invention, 
probe intensities of a sample sequence are compared to statistics about probe intensities of a reference sequence from 
multiple experiments. 

According to another aspect of the invention, a method is disclosed of processing reference and sanple nucleic 
acid sequences to reduce the variations between the experiments by the steps of: 

- providing a plurality of nucleic acid probes; 

• labeling the reference nucleic acid sequence with a first marker; 

- labeling the sample nucleic acid sequence with a second marker; and 
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hybridizing the labeled reference and sample nucleic acid sequences at the same time. 

According to another aspect of the invention, a computer system is used to identify mutations in a sample nucleic 
acid sequence by the steps of: 

5 - inputting a first set of probe intensities, each of the probe intensities in said f isrt set being associated with a nucleic 
acid probe and substantially proportional to the associated nucleic acid probe hybridizing with a reference nucleic 
acid sequence; 

- inputting a second set of probe intensities, each of the probe intensities In said fisrt set being associated with a 
nucleic acid probe and substantially proportional to the associated nucleic acid probe hybridizing with said sample 

10 sequence; 

- the computer system comparing probe intensities in the first set to probe intensities in the second set to 6elect 
hybridization regions where the probe intensities in the first and second sets differ; and 

identifying mutations according to characteristics of the selected regiona 
is According to yet another aspect of the invention, a computer system is used for comparative analysis and visuali- 
zation of multiple sequences by the steps of: 

- displaying at least one reference sequence in a first area on a display device; and 

- displaying at least one sample sequence in a second area on said display device; 

20 

whereby a user is capable of visually comparing the multiple sequences. 

A further understanding of the nature and advantages of the inventions herein may be realized by reference to the 
remaining portions of the specification and the attached drawings. 

25 BRIEF DESCRIPTION OF THE DRAWINGS 

Fig. 1 illustrates an example of a computer system used to execute the software of the present invention; 

Fig. 2 shows a system block diagram of a typical computer system used to execute the software of the present 

invention; 

30 Fig. 3 Illustrates an overall system for forming and analyzing arrays of biological materials such as DNA or RNA; 

Fig. 4 is an illustration of the software for the overall system; 

Rg. 5 illustrates the global layout of a chip formed in the overall system; 

Fig. 6 illustrates conceptually the binding of probes on chips; 

Fig. 7 illustrates probes arranged In lanes on a chip; 
as Fig. 8 illustrates a hybridization pattern of a target on a ch*p with a reference sequence as In Fig. 7; 

Fig. 9 illustrates the high level flow of the intensity ratio method; 

Fig. 1 0A illustrates the high level flow of one implementation of the reference method and Fig. 1 0B shows an analysis 
table for use with the reference method; 

Fig. 1 1 A illustrates the high level flow of another implementation of the reference method; Fig. 11 B shows a data 
40 table for use with the reference method; Fig. 11C 6hows a graph of the normalized sample base intensities minus 
the normalized reference base intensities; and Fig. 11D shows other graphs of data in the data table; 
Fig. 12 illustrates the high level flow of the statistical method; 

Fig. 13 illustrates the pooling processing of a reference and sample nucleic acid sequence; 
Figs. 14A and 14C show graphs of scaled fluorescent intensities of wild-type probes hybridizing with sample and 
45 reference sequences and 14B shows a hypothetical graph of fluorescent intensities of wild-type probes hybridizing 
with two sample sequences and a reference sequence; 

Rg. 1 5 illustrates the high level flow of an embodiment that uses the hybridization data from than one base position 
to identify mutations In a sample sequence; 

Rg. 16 illustrates the main screen and the associated pull down menus for comparative analysis and visualization 
so of multiple experiments; 

Fig. 17 illustrates an Intensity graph window for a selected base; 
Fig. 18 Illustrates multiple intensity graph windows for selected bases; 

Fig. 19 illustrates the intensity ratio method correctly calling a mutation in solutions with varying concentrations; 
Fig. 20 Illustrates the reference method correctly calling a mutant base where the intensity ratio method incorrectly 
68 called the mutant base; and 

Rg. 21 illustrates the output of the ViewSeq™ program with four pretreatment samples and four posttreatment 6am- 
plea 
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DESCRIPTION OF THE PREFERRED EMBODIMENT 
CONTENTS 

6 I. General 

II. Intensity Ratio Method 

III. Reference Method 
IV Statistical Method 
V. Pooling Processing 

10 VI. Comparative Analysis 
VII. Examples 

I. General 

is In the description that follows, the present Invention will be described in reference to a Sun Workstation In a UNIX 
environment The present invention, however, is not limited to any particular hardware or operating system environment 
Instead, those skilled in the art will find that the systems and methods of the present Invention may be advantageously 
applied to a variety of systems, including IBM personal computers running MS-DOS or Microsoft Windows. Therefore, 
the following description of specific systems are for purposes of illustration and not limitation. 

20 Fig. 1 illustrates an example of a computer system used to execute the software of the present Invention. Fig. 1 
shows a computer system 1 which includes a monitor 3, screen 5, cabinet 7, keyboard 9, and mouse 1 1 . Mouse 1 1 may 
have one or more buttons such as mouse buttons 13. Cabinet 7 houses a floppy disk drive 14 and a hard drive (not 
shown) that may be utilized to store and retrieve software programs incorporating the present invention. Although a 
floppy disk 15 is shown as the removable media, other removable tangible media including CD-ROM, flash memory and 

2S tape may be utilized. Cabinet 7 also houses familiar computer components (not shown) such as a processor, memory, 
and the like. 

Fig. 2 shows a system block diagram of computer system 1 used to execute the software of the present invention. 
As in Fig. 1 , computer system 1 includes monitor 3 and keyboard 9. Computer system 1 further includes subsystems 
such as a central processor 52, system memory 54. WO controller 56, display adapter 58, serial port 62, disk 64, network 
so interface 66, and speaker 68. Disk 64 is representative of an Internal hard drive, floppy drive, CD-ROM, flash memory, 
tape, or any other storage medium. Other computer systems suitable for use with the present Invention may include 
additional or fewer subsystems. For example, another computer system could include more than one processor 52 (i.e., 
a multi-processor system) or memory cache. 

Arrows such as 70 represent the system bus architecture of computer system 1 . However, these arrows are ilius- 
35 trative of any interconnection scheme serving to link the subsystems. For example, speaker 68 could be connected to 
the other subsystems through a port or have an Internal direct connection to central processor 52. Computer system 1 
shown in Fig. 2 is but an example of a computer system suitable for use with the present invention. Other configurations 
of subsystems suitable for use with the present Invention will be readily apparent to one of ordinary skill In the art 

The VLSIPS™ technology provides methods of making very large arrays of oligonucleotide probes on very small 
40 chips. See US. Patent No. 5,143,654 and PCT patent publication Nos. WO 90/15070 and 92/10092, each of which is 
incorporated by reference for all purposes. The oligonucleotide probes on the DNA probe array are used to detect 
complementary nucleic acid sequences in a sample nucleic acid of interest (the 'larger nucleic acid). 

The present invention provides methods of analyzing hybridization intensity files for a chip containing hybridized 
nucleic acid probes. In a representative embodiment, the files represent fluorescence data from a biological array, but 
46 the files may also represent other data such as radioactive Intensity data or large molecule detection data. Therefore, 
the present invention is not limited to analyzing fluorescent measurements of hybridizations but may be readily utilized 
to analyze other measurements of hybridization. 

For purposes of illustration, the present invention is described as being part of a computer system that designs a 
chip mask, synthesizes the probes on the chip, labels the nucleic acids, and scans the hybridized nucleic acid probes. 
so Such a system Is fully described in U.S. Patent Application No. 08/249,1 88 which has been Incorporated by reference 
tor all purposes. However, the present invention may be used separately from the overall system for analyzing data 
generated by such systems. 

Fig. 3 illustrates a computerized system tor forming and analyzing arrays of biological materials such as RNA or 
DNA. A computer 1 00 is used to design arrays of biological polymers such as RNA or DNA. The computer 1 00 may be, 
66 tor example, an appropriately programmed Sun Workstation or personal computer or workstation, such as an IBM PC 
equivalent, including appropriate memory and a CPU as shown in Figs. 1 and 2. The computer system 100 obtains 
Inputs from a user regarding characteristics of a gene of interest and other inputs regarding the desired features of the 
array. Optionally, the computer system may obtain information regarding a specific genetic sequence of interest from an 
external or internal database 102 such as GenBank. The output of the computer system 100 is a set of crap design 
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computer files 1 04 in the form of, for example, a switch matrix, as described in PCT application WO 92/10092, and other 
associated computer files. 

The chip design files are provided to a system 106 that designs the lithographic masks used in the fabrication of 
arrays of molecules such as DNA. The system or process 106 may include the hardware necessary to manufacture 
masks 1 1 0 and also the necessary computer hardware and software 1 08 necessary to lay the mask patterns out on the 
mask in an efficient manner. As with the other features in Rg. 3. such equipment may or may not be located at the same 
physical site, but is shown together for ease of illustration in Rg. 3. The system 106 generates masks 110 or other 
synthesis patterns such as chrome-on-glass masks for use in the fabrication of polymer arrays. 

The masks 1 10, as well as selected information relating to the design of the chips from system 100, are used in a 
synthesis system 112. Synthesis system 1 1 2 includes the necessary hardware and software used to fabricate arrays of 
polymers on a substrate or chip 1 14. For example, synthesizer 1 12 includes a light source 1 16 and a chemical flow cell 
1 1 8 on which the substrate or chip 1 14 is placed Mask 1 10 is placed between the light source and the substratefchip. 
and the two are translated relative to each other at appropriate times for deprotection of selected regions of the chip. 
Selected chemical reagents are directed through flow cell 1 1 8 for coupling to deprotected regions* as well as for washing 
and other operations. All operations are preferably directed by an appropriately programmed computer 119, which may 
or may not be the same computer as the computers) used in mask design and mask making. 

The substrates fabricated by synthesis system 1 12 are optionally diced into smaller chips and exposed to marked 
receptors. The receptors may or may not be complementary to one or more of the molecules on the substrata The 
receptors are marked with a label such as a fluorescein label (indicated by an asterisk in Rg. 3) and placed in scanning 
system 120. Scanning system 120 again operates under the direction of an appropriately programmed digital computer 
122, which also may or may not be the same computer as the computers used in synthesis, mask making, and mask 
design. The scanner 120 Includes a detection device 124 such as a confocal microscope or CCD (charge-coupled 
device) that is used to detect the locations where labeled receptor (*) has bound to the substrata The output of scanner 
120 is an image file(s) 124 indicating, in the case of fluorescein labeled receptor, the fluorescence intensity (photon 
counts or other related measurements, such as voltage) as a function of position on the substrata Since higher photon 
counts will be observed where the labeled receptor has bound more strongly to the array of polymers, and since the 
monomer sequence of the polymers on the substrate is known as a function of position, it becomes possible to determine 
the sequence(s) of polymer(s) on the substrate that are complementary to the receptor. 

The image file 124 is provided as input to an analysis system 126 that incorporates the visualization and analysis 
methods of the present Invention. Again, the analysis system may be any one of a wide variety of computer system(s), 
but in a preferred embodiment the analysis system is based on a Sun Workstation or equivalent The present invention 
provides various methods of analyzing the chip design files and the image files, providing appropriate output 12a The 
present invention may further be used to identify specific mutations In a receptor such as DNA or RNA. 

Fig. 4 provides a simplified illustration of the overall software system used in the operation of one embodiment of 
the invention. As shown in Rg. 4, in some cases (such as sequence checking systems) the system first identifies the 
genetic sequenced) or targets that would be of interest in a particular analysis at step 202. The sequences of interest 
may, for example, be normal or mutant portions of a gene, genes that identify heredity, or provide forensic information, 
or be all possible n-mers (where n represents the length of the nucleic acid). Sequence selection may be provided via 
manual input of text files or may be from external sources such as Gen Bank. At step 204 the system evaluates the gene 
to determine or assist the user in determining which probes would be desirable on the chip, and provides an appropriate 
"layout" on the chip for the probes. The chip usually includes probes that are complementary to a reference nucleic acid 
sequence which has a known sequence. A wild-type probe is a probe that will ideally hybridize with the reference 
sequence and thus a wild-type gene (also called the chip wild-type) would ideally hybridize with wild-type probes on the 
chip. The target sequence is substantially similar to the reference sequence except tor the presence of mutations, inser- 
tions, deletions, and the like. The layout implements desired characteristics such as arrangement on the chip that permits 
"reading" of genetic sequence and/or minimization of edge effects, ease of synthesis, and the like. 

Rg. 5 Illustrates the global layout of a chip in a particular embodiment used for sequence checking applications. 
Chip 1 14 Is composed of multiple units where each unit may contain different tilings tor the chip wild-type sequence. 
Unit 1 is shown in greater detail and shows that each unit is composed of multiple ceils which are areas on the chip that 
may contain probes. Conceptually, each unit is composed of multiple sets of related cells. As used herein, the term cell 
refers to a region on a substrate that contains many copies of a molecule or molecules of interest Each unit is composed 
of multiple cells that may be placed in rows (or "lanes") and columns. In one embodiment, a set of five related cells 
includes the following: a wild-type cell 220, "mutation" cells 222, and a "Wank" cell 224. Cefl 220 contains a wild-type 
probe that is the corrplement of a portion of the wild-type sequenca Cells 222 contain "mutation" probes for the wild- 
type sequenca For example, if the wild-type probe Is 3*-ACGT, the probes S'-ACAT, 3'-ACCT, 3'-ACQT, and 3'-AGTT 
may be the "mutation" probes. Cell 224 is the "blank" cell because it contains no probes (also called the "blank" probe). 
As the blank cell contains no probes, labeled receptors should not bind to the chip in this area. Thus, the blank cell 
provides an area that can be used to measure the background intensity. 
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In one embodiment numerous tiling processes are available including sequence tiling, block tiling, and opt-tiling, 
as described below. Of course a wide range of layout strategies may be used according to the Invention herein, without 
departing from the scope of the invention. For example, the probes may be tiled on a substrate in an apparently random 
fashion where a computer system is utilized to keep track of the probe locations and correlate the data obtained from 
e the substrata 

Opt-tiling is the process of tiling additional probes for suspected mutations. As a simple example of opt-tiling, suppose 
the wild-type target sequence is 5 -ACQ WGCA-3* and it is suspected that a mutant sequence has a possible T base 
mutation at the underlined base position. Suppose further that the chip will be synthesized with a "4x3" tiling strategy, 
meaning that probes of four monomers are used and that the monomers in position 3. counting left to right of the probe 
w are varied. 

In opt-tiling, extra probes are tiled for each suspected mutation. The extra probes are tiled as if the mutation base 
is a wild-type base. The following shows the probes that may be generated for this example: 



Table 1 
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Probe Sequences (From 3'- end) 4x3 Opt-Tiling 


Wild 


TGCA 


GCAT 


CATA 


ATAC 


TACG 


A sub. 


TGAA 


GCAT 


CAAA 


ATAC 


TAAG 


Csub. 


TGCA 


GCCT 


CACA 


ATCC 


TACG 


Gsub. 


TGGA 


GCGT 


CAGA 


ATGC 


TAGG 


Tsutx 


TGTA 


GCTT 


CATA 


ATTC 


TATG 


Wild 


TGCA 


GCAA 


CAAA 


AAAC 


AACG 


A sub. 


TGAA 


GCAA 


CAAA 


AAAC 


AAAG 


Csub. 


TGCA 


GCCA 


CACA 


AACC 


AACG 


Gsub. 


TGGA 


GCGA 


CAGA 


AAGC 


AAGG 


Tsub. 


TGTA 


GCTA 


CATA 


AATC 


AATQ 



In the first "chip" above, the top row of the probes (along with one probe below each of the four wild-type probes) should 
bind to the target DNA sequence. However, if the target sequence has a T base mutation as suspected, the labeled 

SB mutant sequence will not bind that strongly to the probes in the columns around column 3. For example, the mutant 
receptor that could bind with the probes in column 2 is 5*-CGTT which may not bind that strongly to any of the probes 
in column 2 because there are T bases at the ends of the receptor and probes (i.e., not complementary). This often 
results in a relatively dark scanned area around a mutation. 

Opt-tiling generates the second "chip" above to handle the suspected mutation as a wild-type base. Thus, the mutant 

40 receptor S'-CGTT should bind strongly to the wild-type probe of column 2 (along with one probe below) and the mutation 
can be further detected. 

Again referring to Fig. 4, at step 206 the masks for the synthesis are designed. At step 208 the software utilizes the 
mask design and layout information to make the DNA or other polymer chips. This software 208 will control relative 
translation of a substrate and the mask, the flow of desired reagents through a flow cell, the synthesis temperature of 

48 the flow cell, and other parameters. At step 21 0, another piece of software is used in scanning a chip thus synthesized 
and exposed to a labeled receptor. The software controls the scanning of the chip, and stores the data thus obtained In 
a file that may later be utilized to extract sequence information. 

At step 212 a computer system according to the present invention utilizes the layout information and the fluorescence 
information to evaluate the hybridized nucleic acid probes on the chip. Among the important pieces of information 

so obtained from probe arrays are the Identification of mutant receptors and determination of genetic sequence of a par- 
ticular receptor. 

Fig. 6 illustrates the binding of a particular target DNA to an array of DNA probes 1 14. As shown In this simple 
example, the following probes are formed in the array (only one probe is shown for the wild-type probe): 

55 
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3 1 -AGAACGT 
AGACCGT 
AGAGCGT 
AGATCGT 



As shown, the set of probes differ by only one base so the probes are designed to determine the Identity of the base at 
that position in the nucleic acid sequence. 

When a fluorescein-labeled (or otherwise marked) target with the sequence 5*-TCTTQCA is exposed to the array, 
ft is complementary only to the probe 3'-AGAACGT, and fluorescein will be primarily found on the surface of the chip 
where 3-AGAACGT is located. Thus, for each set of probes that differ by only one base, the image file will contain four 
fluorescence intensities, one for each proba Each fluorescence intensity can therefore be associated with the base of 
each probe that fe different from the other probes. Additionally, the image file win contain a "blank" cell which can be 
used as the fluorescence intensity of the background. By analyzing the five fluorescence intensities associated with a 
specific base location, it becomes possible to extract sequence information from such arrays using the methods of the 
invention disclosed herein. ^ ^ 

Fig. 7 illustrates probes arranged in lanes on a chip. A reference sequence is shown with five interrogation positions 
marked with number subscripts. An interrogation position is a base position in the reference sequence where the target 
sequence may contain a mutation or otherwise differ from the reference sequence. The cNp may contain five probe cells 
that correspond to each interrogation position. Each probe cell contains a set of probes that have a common base at 
the interrogation position. For example, at the first interrogation position, l lp the reference sequence has a base T. The 
wild-type probe for this interrogation position is 3'-TGAC where the base A in the probe Is complementary to the base 
at the interrogation position in the reference sequence. 

Similarly, there are four "mutant" probe cells tor the first interrogation position. Ii- The four mutant probes are 3- 
TGAC, 3'TGCC, 3'-TGGC, and 3'-TGTC. Each of the four mutant probes vary by a single base at the interrogation 
position. As shown, the wild-type and mutant probes are arranged in lanes on the chip. One of the mutant probes (in 
this case 3'-TGAC) is identical to the wild-type probe and therefore does not evidence a mutation. However, the redun- 
dancy gives a visual indication of mutations as wfll be seen in Fig. 8. 

Still referring to Fig. 7, the chip contains wild-type and mutant probes for each of the other interrogation positions 
l r l B . In each case, the wild-type probe is equivalent to one of the mutant probes. 

Fig. 8 illustrates a hybridization pattern of a target on a chip with a reference sequence as in Fig. 7. The reference 
sequence is shown along the top of the chip tor comparison. The chip includes a WT-lane (wild-type), an A-lane, a C- 
lane, a G-lane, and a T-lane (or U). Each lane is a row of cells containing probes. The cells in the WT-lane contain probes 
that are complementary to the reference sequence. The cells in the A-, C-, G-. and T-lanes contain probes that are 
complementary to the reference sequence except that the named base is at the Interrogation position. 

In one embodiment, the hybridization of probes In a cell is determined by the ftoorescent Intensity (e.g., photon 
counts) of the cell resulting from the binding of marked target sequences. The fluorescent intensity may vary greatly 
among cells. For simplicity, Fig. 8 shows a high degree of hybridization by a cell containing a darkened area. The WT- 
lane allows a simple visual Indication that there Is a mutation at interrogation position U because the wiJd-type cell Is not 
dark at that position. The cell in the C-iane Is darkened which Indicates that the mutation is from T->G (mutant probe 
ceils are complementary so the C-cell indicates a G mutation). 

In practice, the fluorescent intensities of ceils near an interrogation position having a mutation are relatively dark 
creating "dark regions" around a mutation. The lower fluorescent Intensities result because the cells at Interrogation 
positions near a mutation do not contain probes that are perfectly complementary to the target sequence; thus, the 
hybridization of these probes with the target sequence is lower. For example, the relative Intensity of the cells at Inter- 
rogation positions l 3 and l 6 may be relatively tow because none of the probes therein are complementary to the target 
sequence. 
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For ease of reference, one may call bases by assigning the bases the following codes: 



Code 


Group 


Meaning 


A 


A 


Adenine 


C 


C 


Cytosine 


Q 


Q 


Guanine 


T 


T(U) 


Thymine (Uracil) 


M 


AorC 


aMtno 


R 


AorG 


puRJne 


W 


AorT(U) 


Weak interaction (2 H bonds) 


Y 


CorT(U) 


pYrimidine 


S 


C or G 


Strong interaction (3 H bonds) 


K 


GorT(U) 


Keto 


V 


A, CorG 


notT(U) 


H 


A.CorT(U) 


notG 


D 


A,GorT<U) 


notC 


B 


CGorTXU) 


not A 


N 


A,C.G,orT(U) 


insufficient intensity to call 


X 


A,C,G,orT(U) 


Insufficient discrimination to call 



Most of the codes conform to the lUPAC standard. However, code N has been redefined and code X has been added. 
H. Intensity Ratio Method 

The Intensity ratio method is a method of calling bases in a sample nucleic acid sequence. The intensity ratio method 
is most accurate when there is good discrimination between the fluorescence intensities of hybrid matches and hybrid 
mismatches. H there is insufficient discrimination, the intensity ratio method assigns a corresponding ambiguity code to 
the unknown base. 

For simplicity, the intensity ratio method will be described as being used to identify one unknown base in a sample 
nucleic acid sequence. In practice, the method is used to identify many or all the bases in a nudelc acid sequence. 

The unknown base will be identified by evaluation of up to four mutation probes and a 'blank" cell, which is a location 
where a labeled receptor should not bind to the chip since no probe is present. For example, suppose a DNA sequence 
of interest or target sequence contains the sequence 5-AGAA&CTGC-3' with a possible mutation at the underlined base 
positioa Suppose that 5-mer probes are to be synthesized for the target sequence. A representative wild-type probe of 
S'-TTGGA Is complementary to the region of the sequence around the possible mutation. The "mutation" probes will be 
the same as the wild-type probe except for a different base at the third position as follows: 3-TTAGA, 3'-TTCQA, 3 « 
TTGGA, and 3'-TTTGA. 

If the fluoresce rrtiy marked sample sequence Is exposed to the above four mutation probes, the Intensity should be 
highest for the probe that binds most strongly to the sample sequence Therefore, if the probe 3-TTTGA shows the 
highest intensity, the unknown base In the sample win generally be called an A mutation because the probes are com- 
plementary to the sample sequence. 

The mutation probes are identical to the wild-type probes except that they each contain one of the four A, 0,0, or 
T "mutations" for the unknown base. Although one of the "mutation" probes will be identical to the wild-type probe, such 
redundant probes are Intentionally synthesized for quality control and design consistency. 

The identity of the unknown base Is preferably determined by evaluating the relative fluorescence Intensities of up 
to four of the mutation probes, and the "blank" ceil. Because each mutation probe is identifiable by the mutation base, 
_a mutation probe's intensity will be referred to qs the T^ase intensity" of the mutation base. 

As a simple example of the intensity ratio method, suppose a gene of interest (target) is an HIV protease gene with 
the sequence 5 , -ATGTGGAGAGTTGTA-3' (SEQ ID NO:1). Suppose further that a sample sequence is suspected to 
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have the same sequence as the target sequence except for a mutation of base C to base T at the underlined base 
position. Although hundreds of probes may be synthesized on the chip, the complementary mutation probes synthesized 
to detect a mutation in the sample sequence at the suspected mutation position may be as follows: 

3-TATC 

3-TCTC 

3'-TGTC (wild-type) 
3'-TTTC 

The mutation probe 3'-TGTC is also the wild-type probe as it should bind most strongly with the target sequence. 

After the sample sequence is labeled, hybridized on the chip, and scanned, suppose the following fluorescence 
intensities were obtained: 

3'-TATC->45 

3'-TCTC->8 

3'-TGTC ->32 

3'-TTTC -> 12 

where the intensity is measured by the photon count detected by the scanner. The "Wank" cell had a fluorescence 
intensity of 2. The photon counts in the examples herein are representative (not actual data) and provided for illustration 
purposes. In practice, the actualphoton counts will vary greatly depending on the experiment parameters and the scanner 
utilized. 

Although each fluorescence intensity is from a probe, the probes may be characterized by their unique mutation 
base so the bases may be said to have the following intensities: 
A->45 
C->8 
G->32 
T->12 

Thus, base A will be described as having an intensity of 45, which corresponds to the intensity of the mutation probe 
with the mutation base A. 

Initially, each mutation base intensity is reduced by the background or Hank" cell intensity. This is done as follows: 

A-> 45-2-43 
C->8-2 = 6 
Q-> 32-2-30 
T-> 12-2-10 

Then, the base Intensities are sorted in descending order of intensity. The above bases would be sorted as follows: 
A->43 
G->30 
T->10 
C->6 

Next the highest intensity base is compared to the second highest intensity base. Thus, the ratio of the intensity of base 
A to the intensity of base G is calculated as follows: A.*G » 43 / 30 - 1 .4. The ratio A:G is then compared to a predetermined 
ratio cutoff, which is a number that specifies the ratio required to identify the unknown base. For example, if the ratio 
cutoff is 1.2, the ratio AG is greater than the ratio cutoff (1.4 > 1.2) and the unknown base Is called by the mutation 
probe containing the mutation A. As probes are complementary to the sample sequence, the sample sequence is called 
as having a mutation T, resulting In a called sample sequence of 5WGTGGAI*GTTGTA-3 , (SEQ ID NO:2). 

As another example, suppose everything else is the same as in the previous example except that the sorted back- 
ground adjusted intensities were as follows: 

C->42 

A->40 

G-> 10 

T->8 

The ratio of the highest intensity base to the second highest intensity base (C:A) Is 1 .OS Because this ratio is not greater 
than the ratio cutoff of 1 .2, the unknown base will be called as being ambiguously one of two or more bases as follows. 

The second highest intensity base is then compared to the third highest base. The ratio of A:G is 4. The ratio of A:G 
is then compared to the ratio cutoff of 1 .2. As the ratio A:G is greater than the ratio cutoff (4 > 1 2), the unknown base 
is caned by the mutation probes containing the mutations C or A. As probes are complementary to the sample sequence, 
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the sample sequence is called as having either a mutation Q or T, resulting in a sample sequence of 5'-ATGTGQA& 
AGTTGTA-3 1 (SEQ ID N03) where K is the IUPAC code for G or T(U). 

The ratio cutoff in the previous examples was equal to 1 .2, However, the ratio cutoff will generally need to be adjusted 
to produce optimal results for the specific chip design and wild-type target Also, although the ratio cutoff used has been 
the same for each ratio comparison, the ratio cutoff may vary depending on whether the ratio comparisons involve the 
highest second highest third highest etc. intensity base. 

Fig. 9 illustrates the high level flow of the intensity ratio method At step 302 the four base intensities are adjusted 
by subtracting the background or "blank" cell intensify from each base intensity. Preferably, if a base intensity is then 
less than or equal to zero, the base intensity is set equal to a small positive number to prevent division by zero or negative 
numbers in future calculations. 

At step 304 the base intensities are sorted by intensity. Each base is then associated with a number from 1 to 4. 
The base with the highest intensity is 1 . second highest 2, third highest 3, and fourth highest 4. Thus, the Intensity of 
base 1 £ base 2 & base 3 & base 4. 

At step 306 the highest intensity base (base 1) is checked to see if it has sufficient intensity to call the unknown 
base. The intensity is checked by determining rf the intensity of base 1 is greater than a predetermined background 
difference cutoff. The background difference cutoff is a number that specifies the intensity a base intensity must be over 
the background intensity in order to correctly call the unknown base. Thus, the background adjusted base intensity must 
be greater than the background difference cutoff or the unknown is not callable. 

If the intensity of base 1 is not greater than the background difference cutoff, the unknown base is assigned the 
code N (insufficient intensity) as shown at step 308. Otherwise, the ratio of the intensity of base 1 to base 2 is calculated 
as shown at step 310. 

At step 312 the ratio of intensity of bases 1 ;2 is compared to the ratio cutoff. If the ratio 1 :2 is greater than the ratio 
cutoff, the unknown base is called as the complement of the highest Intensity base (base 1) as shown at step 314. 
Otherwise, the ratio of the intensity of base 2 to base 3 is calculated as 6hown at step 316. 

At step 31 8 the ratio of intensity of bases 2:3 is compared to the ratio cutoff. H the ratio 2:3 Is greater than the ratio 
cutoff, the unknown base is called as being ah ambiguity code specifying the complements of the highest or second 
highest intensity bases (base 1 or 2) as shown at step 320. Otherwise, the ratio of the intensity of base 3 to base 4 is 
calculated as shown at step 322. 

At step 324 the ratio of intensity of bases 3:4 is compared to the ratio cutoff. If the ratio 3:4 Is greater than the ratio 
cutoff, the unknown base is called as being an ambiguity code specifying the complements of the highest second 
highest or third highest bases (base 1. 2 or 3) as shown at step 326. Otherwise, the unknown base is assigned the 
code X (insufficient discrimination) as shown at step 328. 

The advantage of the intensity ratio method is that it Is very accurate when there Is good discrimination between 
the fluorescence intensities of hybrid matches and hybrid mismatches. However, if the base corresponding to a correct 
hybrid gives a lower intensity than a mismatch (e.g„ as a result of cross-hybridization), incorrect identification of the 
base will resuK. For this reason, however, the method is useful for comparative assessment of hybridization quality and 
as an indicator of sequence-specific problem spots. For example, the intensity ratio method has been used to determine 
that ambiguities and miscalls tend to be very different from sequence to sequence, and reflect predominantly the com- 
position and repetitiveness of the sequence. It has also been used to assess improvements obtained by varying hybrid- 
ization conditions, sample preparation, and post-hybridization treatments (e.g.. RNase treatment). 

III. Reference Method 

The reference method is a method of calling bases in a sample nucleic acid sequence. The reference method 
depends very little on discrimination between the fluorescence intensities of hybrid matches and hybrid mismatches, 
and therefore is much less sensitive to cross-hybridization. The method compares the probe intensities of a reference 
sequence to the probe intensities of a sample sequence. Any significant changes are flagged as possible mutations. 
There are two implementations of the reference method disclosed herein. 

For simplicity, the reference method will be described as being used to identify one unknown base in a sample 
nucleic acid sequence. In practice, the method is used to Identify many or all the bases in a nucleic acid sequence. 

The unknown base will be called by comparing the probe intensities of a reference sequence to the probe Intensities 
of a sample sequenca Preferably, the probe intensities of the reference sequence and the sample sequence are from 
chips having the same chip wild-type. However, the reference sequence may or may not be exactly the same as the chip 
wild-type, as it may have mutations. 

The bases at the same position in the reference and sample sequences will each be associated with up to four 
mutation probes and a "blank" cell. The unknown base In the sample sequence is called by comparing probe intensities 
of the sample sequence to probe intensities of the reference sequence. For example, suppose the chip wild-type contains 
the sequence 5-AQACCTTGC-3' and it is suspected that the sample has a poseHe mutation at the underlined base 
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position, which is the unknown base that will be called by the reference method. The •mutation" probes for the sample 
sequence may be as follows: 3'-GAAA, 3'-GCAA, 3*-GGAA, and 3'-GTAA, where 3'-GGAA is the wild-type probe. 

Suppose further that a reference sequence, which differs from the chip wild-type by one base mutation, has the 
sequence S-AGACATTGC-tf where the mutation base is underlined. The "mutation" probes for the reference sequence 

s may be as follows: 3**TGAAA, 3'-TGCAA, 3'-TGGAA. and 3'-TGTAA, where 3'-TGTAA is the reference wild-type probe 
since the reference sequence is known. Although generally the sample and reference sequences were tiled with the 
same chip wild-type, this is not required, and the tiling methods do not have to be identical as shown by the use of two 
probe lengths in the example. Thus, the unknown base will be called by comparing the "mutation" probes of the sample 
sequence to the "mutation" probes of the reference sequence. As before, because each mutation probe is identifiable 

io by the mutation base, the mutation probes' intensities will be referred to as me "base intensities" of their respective 
mutation base6. 

As a simple example of one implementation of the reference method, suppose a gene of interest (target) has the 
sequence 5'-AAAACTGAAAA-3* (SEQ ID NO:4). Suppose a reference sequence has the sequence S-AAAACGGAAAA- 
3* (SEQ ID NO:5). which differs from the target sequence by the underlined base. The reference sequence is marked 
15 and exposed to probes on a chip with the target sequence being the chip wild-type. Suppose further that a sample 
sequence is suspected to have the same sequence as the target sequence except for a mutation at the underlined base 
position in 5-AAAACIGAAAA-3* (SEQ ID NO:4). The sample sequence is also marked and exposed to probes on a 
chip with the target sequence being the chip wild-type. After hybridization and scanning, the following probe intensities 
(not actual data) were found for the respective complementary probes: 

20 



25 



Reference 


Sample 


3*-TGAC->12 
3'-TGCC ->9 
3'-TGGC ->80 
3'-TGTC->15 


3'-GACT->11 
3M3CCT->30 
3'-GGCT ->60 
3'-GTCT->6 



30 

Although each fluorescence intensity is from a probe, the probes may be identified by their unique mutation base so the 
bases may be said to have the following intensities: 



Reference 


Sample 


A->12 


A->11 


C->9 


C->30 


G->80 


G->60 


T->15 


T->6 



Thus, base A of the reference sequence will be described as having an intensity of 1 2, which corresponds to the intensity 
of the mutation probe with the mutation base A. The reference method win now be described as calling the unknown 
base in the sample sequence by using these intensities. 

Fig. 1 0A illustrates the high level flow of one implementation of the reference method. For illustration purposes, the 

so reference method is described as filling in the columns (identified by the numbers along the bottom) of the analysis table 
shown In Fig. 10B. However, the generation of an analysis table Is not necessary to practice the method. The analysis 
table is shown to aid the reader in understanding the method. 

At step 402 the four base intensities of the reference and sample sequences are acljusted by subtracting the back- 
ground or "blanl^ cell intensity from each base intensity. Each set of "mutation" probes has an associated "blankT cell. 

55 Suppose that the reference "Wank" ceil intensity Is 1 and the sample "Wank" cell intensity is 2. The base intensities are 
then background subtracted as follows: 
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Reference 


Sample 


A->12-1«H 
C -> 9 - 1 - 8 
G-> 80-1-79 
T->15-1«14 


A->11 -2-9 
C-> 30-2-28 
G-> 60-2-58 
T->6-2-4 



Preferably, If a base intensity is then less than or equal to zero, the base intensity is set equal to a small positive nuntoer 
to prevent division by zero or negative numbers in future calculations. 

For identification, the position of each base of interest in the reference and sample sequences is placed in colurm 
1 of the analysis table. Also, since the reference sequence is a known sequence, the base at this position is known and 
is referred to as the reference wild-type. The reference wild-type Is placed in column 2 of the analysis table, which Is C 
for this example. 

At step 404 the base intensity associated with the reference wild-type (column 2 of the analysis table) is checked 
to see « It has sufficient intensity to call the unknown base. In this example, the reference wild-type is C. However, the 
base intensity associated with the wild-type is the Q base intensity, which is 79 in this example. This is because the base 
intensities actually represent the complementary "mutation" probes. The Q base intensity Is checked by determining If 
its intensity is greater than a predetermined background difference cutoff. The background difference cutoff is a number 
that specifies the intensity the base intensities must be above the background intensity in order to correctly call the 
unknown base. Thus, the base intensity associated with the reference wild-type must be greater than the background 
difference cutoff or the unknown base is not callable. 

If the background difference cutoff is 5, the base intensity associated with the reference wild-type has sufficient 
intensity (79 > 5) so a P (pass) is placed in column 3 of the analysis table as shown at step 406. Otherwise, at step 407 
an F (fail) is placed in column 3 of the analysis table. 

At step 408 the ratio of the base intensity associated with the reference wild-type to each of the possible bases are 
calculated. The ratio of the base intensity associated with the reference wild-type to itself will be 1 and the other ratios 
will usually be greater than 1 . The base intensity associated with the reference wild-type is Q so the following ratios are 
calculated: 

G:A-> 79/11 -7.2 
Q:C -> 79/8 -9.9 
G:G-> 79/79 -1.0 
G:T-> 79/14- 5.6 

These ratios are placed in columns 4 through 7 of the analysis table, respectively 

At step 410 the highest base intensity associated with the sample sequence is checked to see if it has sufficient 
intensity to call the unknown base. The highest base intensity is checked by determining If the Intensity Is greater than 
the background difference cutoff. Thus, the highest base Intensity must be greater than the background difference cutoff 
or the unknown base is not callable. 

Again, H the background difference cutoff is 5, the highest base Intensity, which is G in this example, has sufficient 
intensity (58 > 5) so a P (pass) is placed in column 8 of the analysis table as shown at step 41 2. Otherwise, at step 413 
an F (tali) is placed in column 8 of the analysis table. 

At step 414 the ratios of the highest base intensity of the sample to each of the possible bases are calculated. The 
ratio of the highest base intensity to itself will be 1 and the other ratios wfll usually be greater than 1. Thus, the highest 
base intensity is Q so the following ratios are calculated: 

G:A-> 58/9-6.4 

G:C-> 58/28 - 2.3 

G:G-> 58/58 -1.0 
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Q:T-> 58/4 o 14.5 

These ratios are placed in columns 9 through 12 of the analysis table, respectively. 

At step 416 H both the reference and sample sequence probes failed to have sufficient intensity to call the unknown 
base, meaning there is an F in columns 3 and 8 of the analysis table, the unknown base is assigned the code N 
(insufficient intensity) as shown at step 418. An 'N' is placed in column 17 of the analysis table. Additionally, a confidence 
code of 9 is placed in column 18 of the analysis table where the confidence codes have the following meanings: 



Code 
"~0 

1 

2 

3 

4 

5-8 
9 



Meaning 

Probable reference wild-type 
Probable mutation 

Reference sufficient intensity, insufficient intensity in sample suggests possible mutation 
Borderline differences, unknown base ambiguous 

Sample sufficient intensity, insufficient intensity in reference to allow comparison 
Currently unassigned 

Insufficient intensity in reference and sample, no interpretation possible 



The confidence codes are useful tor indicating to the user the resulting analysis of the reference method. 

Atstep 420 if only the reference sequence probes failed to have sufficient intensity to callthe unknown base, meaning 
there is an'F in column3 and a F in column 8 of the analysis table, the unknown base is assigned the code N (insufficient 
intensity) as shown at step 422. An 'N' is placed in column 17 and a confidence code of 4 is placed in column 18 of the 
analysis table. 

At step 424 if only the sample sequence probes failed to have sufficient intensity to call the unknown base, meaning 
there Is a F in column 3 and a F In column 8 of the analysis table, the unknown base is assigned the code N (insufficient 
intensity) as shown at step 426. An 'N* is placed in column 17 and a confidence code of 2 is placed in column 18 of the 
analysis table. 

In this example, both the reference and sample sequence probes have sufficient intensity to call the unknown base. 
At step 428 the ratios of the reference ratios to the sample ratios for each base type are calculated. Thus, the ratio AA 
(column 4 to column 9) is placed in column 13 of the analysis table. The ratio C:C (column 5 to column 10) is placed in 
column 14 of the analysis table. The ratio G:G (column 6 to column 1 1) is placed in column 15 of the analysis table. 
Lastly, the ratio T:T (column 7 to column 12) is placed in column 16 of the analysis table. These ratios are calculated as 
follows: 

A:A->7.2/6.4-1.1 
C:C -> 9.9/2.3-4.3 
G:G-> 1.0/1.0- 1.0 
T:T-> 5.6/14.5-0.4 

The unknown base is called by comparing these ratios of ratios to two predetermined values as follows. 

At step 430 if all the ratios of ratios (columns 13 to 16 of the analysis table) are less than a predetermined lower 
ratio cutoff, the unknown base is assigned the code of the reference wild-type as shown at step 432. Thus, the code tor 
the reference wild-type (as shown in column 2) would be placed in column 17 and a confidence code of 0 would be 
placed in column 18 of the analysis table. . 

At step 434 if aO the ratios of ratios are less than a predetermined upper ratio cutoff, the unknown base is assigned 
an ambiguity code that indicates the unknown base may be any one of the bases that has a complementary ratio of 
ratios greater than the lower ratio cutoff and less than the upper ratio cutoff as shown at step 436. Thus, if the ratio of 
ratios for A A C:C and G:G are all greater than the lower ratio cutoff and less than the upper ratio cutoff, the unknown 
base would be assigned the code B (meaning "not A"). This is because the ratios of ratios are complementary to their 
respective base as follows: 
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A:A->T 
C:C->Q 
G:G->C 

so the unknown base would be called as being either C. G, or T, which Is identified by the IUPAC code B. This ambiguity 
code would be placed in column 17 and a confidence code of 3 would be placed incolumn 18 of the analysis tabla 

At step 438 at least one of the ratios of ratios Is greater than the upper ratio cutoff and the unknown base is caned 
as the base complementary to the highest ratio of ratios. The code for the base complementary to the highest ratio of 
ratios would be placed in column 17 and a confidence code of 1 would be placed in column 18 of the analysis tabla 

Assume tor the purposes of this example that the lower ratio cutoff to 1 .5 and the upper ratio cutoff is 3. Again, the 
ratios of ratios are as follows: 

A*->1.1 
C:C->4.3 
G:G ->1.0 
T:T-> 0.4 

As all the ratios of ratios are not less than the upper ratio cutoff, the unknown base is called the base complementary 
to the highest ratio of ratios. The highest ratio of ratios is C:C. which has a complementary base Q. Thus, the unknown 
base Is called Q which is placed in column 1 7 and a confidence code of 1 is placed In column f 8 of the analysis tabla 

The example shows hew the unknown base in the sample nucleic acid sequence was correctly called as base Q. 
Although the complementary "mutation" probe associated with the base Q (3'-GCCT) did not have the highest fluores- 
cence intensity, the unknown base was called as base G because the associated "mutation" probe had the highest ratio 
increase over the other "mutation" probea 

Fig. 11 A illustrates the high level flow of another implementation of the reference method. As in the previous Imple- 
mentation, this implementation also compares the probe intensities of a reference sequence to the probe intensities of 
a sample sequence. However, this implementation differs conceptually from the previous Implementation in that neigh* 
boring probe intensities are also analyzed, resulting in more accurate base caning. 

As a simple example of this implementation of the reference method, suppose a reference sequence has a sequence 
of S'-AAACCCAATCCACATCA-S' (SEQ ID NO:6) and a sample sequence has a sequence of 5'-AAACCCAQTCCA- 
CATCA-3' (SEQ ID NOT), where the mutant base is underlined. Thus, there is a mutation of A to G. Suppose further 
that the reference and sample sequences are tiled on chips with the reference sequence being the chip wild-typa This 
implementation of the reference method will be described as identifying this mutation basa 

For illustration purposes, this implementation of the reference method Is described as filling in a data table shown 
in Fig. 1 1 B (SEQ ID NO:6, SEQ ID NO:28, SEQ ID NO:29). Although the data table contains more data than is required 
for this implementation, the portions of the data table that are produced by steps in Rg. 1 1 A are shown with the same 
reference numerals. "Die generation of a data table is not necessary, however, and is shown to aid the reader In under- 
standing the method. The mutant base position is at position 241 In the reference and sample sequences, which is 
shown in bold in the data tabla 

At step 502 the base intensities of the reference and sample sequences are adjusted by subtracting the background 
or "blank" ceil intensity from each base intensity. Preferably, if a base intensity is then less than or equal to zero, the 
base Intensity is set equal to a small positive number to prevent division by zero or negative numbers. In the data table, 
data 502A is the background subtracted base intensities for the reference sequence and data 502B is the background 
subtracted base intensities for the sample sequence (also called the "mutant" sequence In the data table). 

At step 504 the base Intensity associated with the reference wild-type Is checked to see if It has sufficient Intensity 
to call the unknown basa In this example, the reference wild-type Is base A at position 241 . The base Intensity associated 
with the reference wild-type is Identified by a lower case "a" In the left hand column. Thus, the base Intensities in the 
data table are not identified by their complements and the reference wild-type at the mutation position has an intensity 
of 385. The reference wild-type intensity of 385 is checked by determining if its intensity Is greater than a predetermined 
background difference cutoff. The background difference cutoff Is a number that specifies the Intensity the base Intensities 
must be over the background intensity in order to correctly call the unknown basa Thus, the base intensity associated 
-with the reference wild-type-must be greater than the background difference cutoff or the unknown base is not cailaHa 
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ff the base intensity associated with the reference wild-type is not greater than the background difference cutoff, the 
wild-type sequence would fail to have sufficient intensity as shown at step 506. Otherwise, at step 508 the wild-type 
sequence would pass by having sufficient intensity. 

At step 510 calculations are performed on the background subtracted base intensities of the reference sequence 
in order to "normalize* the intensities. Each position in the reference sequence has tour background subtracted base 
intensities associated with it The ratio of the intensity of each base to the sum of the intensities of the possible bases 
(all four) is calculated, resulting in four ratios, one for each base as shown in the data table. Thus, the following ratios 
would be calculated at each position in the reference sequence: 

A ratio- A/(A + C + G + T) 

Cratio»C/(A + C + Q + T) 

G ratio- G/(A + C+G + T) 

T ratio »T/(A + C + G + T) 

At position 241, A ratio would be the wild-type ratia These ratios are generally calculated In order to "normalize" the 
intensity data as the photon counts may vary widely from experiment to experiment Thus, the ratios provide a way of 
reconciling the Intensity variations across experiments. Preferably, if the photon counts do not vary widely from experi- 
ment to experiment, the probe intensities do not need to be "normalized." 

At step 512 the highest base intensity associated with the sample sequence Is checked to see If it has sufficient 
intensity to call the unknown base. The intensity is checked by determining if the highest intensity sample base Is greater 
than the background difference cutoff. If the intensity is not greater than the background difference cutoff, the sample 
sequence fails to have sufficient Intensity as shown at step 514. Otherwise, at step 516 the sample sequence passes 
by having sufficient intensity. 

At step 518 calculations are performed on the background subtracted base intensities of the sample sequence in 
order to "normalize" the intensities. Each position in the sample sequence has four background subtracted base inten- 
sities associated with it. The ratios of the intensity of each base to the sum of the intensities of the possible bases (all 
four) are calculated, resulting In four ratios, one for each base as shown in the data table. 

At step 520 H either the reference or sample sequences failed to have sufficient intensity, the unknown base is 
assigned the code N (insufficient intensity) as shown at step 522. 

At step 524 the normalized base intensities of the reference sequence are subtracted from the normalized base 
intensities of the sample sequence. Thus, at each position the following calculations are performed: 

A Difference • Sample A Ratio - Reference A Ratio 

C Difference - Sample C Ratio - Reference C Ratio 

G Difference - Sample G Ratio - Reference Q Ratio 

T Difference = Sample T Ratio - Reference T Ratio 

where the reference and sample ratios are calculated at steps 510 and 518, respectively. The base differences resulting 
from these calculations are shown in the data table. 

At step 526 each position Is checked to see if there is a base difference greater than an upper difference cutoff and 
a base difference lower than a lower difference cutoff. For example. Fig. 11C shows a graph the normalized sample 
base intensities minus the normalized reference base Intensities. Suppose that the upper difference cutoff Is 0.1 5 and 
the lower difference cutoff is -0.15 as shown by the horizontal lines in Fig. 1 1C. At the mutation position (labded with a 
reference 0), the G difference is 0.28 which is greater than 0.15, the upper difference cutoff. Similarly, the A difference 
is -0.32 which is less than -0.15. the lower difference cutoff. As there is a base difference above the upper difference 
cutoff and a base difference below the lower difference cutoff, there may be mutation at this position. 

If there is neither a base difference above the upper difference cutoff nor a base difference below the lower difference 
cutoff, the base at that position is assigned the code of the reference wild-type base as shown at step 52a 

At step 530 the ratio ol the highest background subtracted base intensity in the sample to the background subtracted 
reference wild-type base intensity is calculated. For example, at the mutation position 241 in the data table, the highest 
background subtracted base intensity in the sample is 571 (base G). The background subtracted reference wild-type 
base intensity is 385 (base A). The ratio of 571 :385 is calculated and results in 1 .48 as shown in the data table. 
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At step 532 these ratios are compared to a ratio at a neighboring position. The ratio for the n m position is subtracted 
from the ratio for the r* position, where r ■ n + 1 . For example, at the mutation position 241 in the data table, the ratio 
at position 242 (which equals 1 .02) is subtracted from the ratio at position 241 (which equals 1 .48). It has been found 
that a mutant can be confidently detected by analyzing the difference of these neighboring ratios. 

Fig. 1 1 D shows other graphs of data In the data table. Of particular importance Is the graph identified as 532 because 
this is a graph of the calculations at step 532. The pattern 6hown in a box in graph 532 has been found to be characteristic 
of a mutation. Thus, if this pattern is detected, the base is called as the base (or bases) with a normalized difference 
greater than the upper difference cutoff as shown at step 536. For example, the pattern was detected and at step 526 
it was shown that base Q had a normalized difference of 0.28, which is greater than the upper difference cutoff of 0.15. 
Therefore, the base at position 241 in the sample sequence is called a base G, which is a mutation from the reference 
sequence (A to Q). 

If the pattern is not detected at step 534, the base at that position is assigned the code of the reference wild-type 
base as shewn at step 538. 

This second implementation of the reference method is preferable in some instances as it takes Into account probe 
intensities of neighboring probes. Thus, the first implementation may not have detected the A to Q mutation in this 
example. - . . 

The advantage of the reference method is that the correct base can be called even in the presence of significant 
levels of cross-hybridization, as long as ratios of intensities are fairly consistent from experiment to experiment. In prac- 
tice, the number of miscalls and ambiguities is significantly reduced, while the number of correct calls is actually 
increased, making the reference method very useful for identifying candidate mutations. The reference method has also 
been used to compare the reproducibility of experiments in terms of base calling. 

iv. Statistical Method 

The statistical method is a method of calling bases in a sample nucleic acid sequence The statistical method utiOzes 
the statistical variation across experiments to call the bases. Therefore, the statistical method is preferable when data 
from multiple experiments is available and the data is fairly consistent across the experiments. The method compares 
the probe intensities of a sample sequence to statistics of probe intensities of a reference sequence in rnuffiple experi- 
ments. 

For simplicity, the statistical method will be described as being used to identify one unknown base in a sample 
nucleic acid sequence In practice, the method is used to identify many or all the bases In a nucleic acid sequence. 

The unknown base wfll be called by comparing the probe intensities of a sample sequence to statistics on probe 
intensities of a reference sequence in multiple experiments. Generally, the probe intensities of the sample sequence 
and the reference sequence experiments are from chips having the same chip wild-type. However, the reference 
sequence may or may not be equal to the chip wild-type, as It may have mutations. 

A base at the same position in the reference and sample sequences will be associated with up to four mutation 
probes and a tlank" cell. As before, because each mutation probe is identifiable by the mutation base, the mutation 
probes* intensities win be referred to as the "base intensities" of their respective mutation bases. 

As a simple example of the statistical method, suppose a gene of interest (target) has the sequence 5*-AAAACT- 
GAAAA-3' (SEQ ID NO:4). Suppose a reference sequence has the sequence 5'-AAAAC£.GAAAA-3' (SEQ ID NO;5), 
which differs from the target sequence by the underlined base. Suppose further that a sample sequence is suspected 
to have the same sequence as the target sequence except for a T base mutation at the underlined base position In 5'- 
AAAACIGAAAA-3' (SEQ ID NO*). Suppose that in multiple experiments the reference sequence Is marked and 
exposed to probes on a chip. Suppose further the sample sequence is also marked and exposed to probes on a chfe. 

The following are complementary -mutation" probes that could be used tor a reference experiment and the sample 
sequence: 



Reference 


Sample 


3-TGAC 


3'-GACT 


3'-TGCC 


3--GCCT 


ff-TOGC 


y-GQCT 


3'-TGTC 


3*-QTCT 



The "mutation" probes shown for the reference sequence may be from only one experiment the other experiments may 
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have different Watton" probes, chip wild-types, tiling methods, and the like. Although each fluorescence intensity is 
from a probe, since the probes may be identified by their unique mutation bases, the probe intensities may be identified 
by their respective bases as follows: 



Reference 


Sample 


3'-TGAC -> A 
3'-TGCC -> C 
3*-TGGC ->G 
3VTGTC->T 


3'-GACT -> A 
3'-GCCT->C 
3M3GCT->G 
3'-GTCT->T 



Thus base A of the reference sequence will be described as having an intensity which corresponds to the Intensity of 
the mutation probe with the mutation base A. The statistical method will now be described as calling the unknown base 
in the sample sequence by using this example. 

Fig. 12 illustrates the high level flow of the statistical method. At step 602 the four base intensities associated with 
the sample sequence and each of the multiple reference experiments are adjusted by subtracting the background or 
"blank" cell Intensity from each base intensity. Preferably, H a base intensity is then less than or equal to zero, the base 
intensity is set equal to a small positive number to prevent division by zero or negative numbers. 

At step 604 the intensities of the reference wild-type bases in the multiple experiments are checked to see if they 
all have sufficient intensity to call the unknown basa The intensities are checked by determining if the intensity of the 
reference wild-type base of an experiment is greater than a predetermined background difference cutoff. The wild-type 
probe shown earlier tor the reference sequence is 3'-TGGC. and trnjs the G base irtensity is the wild-type base intense 
These steps are analogous to steps in the other two methods described herein. 

If the intensity of any one of the reference wild-type bases is not greater than the background difference cutoff, the 
wild-type experiments fail to have sufficient intensity as shown at step 606. Otherwise, at step 608 the wild-type exper- 
iments pass by having sufficient Intensity. 

At step 610 calculations are performed on the background subtracted base intensities of each of the reference 
experiments in order to "normalize" the intensities. Each reference experiment has four background subtracted base 
intensities associated with H: one wild-type and three tor the other possible bases, in this example, the G base Intensity 
is the wild-type, the A, C, and T base intensities being the "other" intensities. The ratios of the intensity of each base to 
the sum of the intensities of the possfole bases (all four) are calculated, giving one wild-type ratio and three "other" ratios. 
Thus, the following ratios would be calculated: 

Aratio-A/(A + C + G + T) 
Cratio-C/(A + C + G + T) 
Gratio- G/(A + C + G + T) 



T ratio «T/(A + C+Q + T) 

where G ratio is the wild-type ratio and A, C. and T ratios are the "other" ratios. These four ratios are calculated tor each 
reference experiment. Thus if the number of reference experiments is n, there would be 4n ratios calculated. These 
ratios are generally calculated in order to "normalize" the intensity data, as the photon counts may vary widely from 
experiment to experiment However. H the probe intensities do not vary widely from experiment to experiment the probe 

intensities do not need to be "normalized." 

At s tep 61 2 statistics are prepared for the ratios calculated for each of the reference experiments, As stated before, 
each reference experiment will be associated with one wild-type ratio and three "other ratios. The mean and standard 
deviation are calculated tor all the wild-type ratios. The mean and standard deviation are also calculated tor each of the 
other ratios, resulting in three other means and standard deviations tor each of the bases that is not the wild-type base. 
Therefore, the following would be calculated: 



Meanand standard deviation of Arafioe 



Mean and standard deviation of C ratios 



17 



EP 0717 113 A2 



Mean and standard deviation of Q ratios 



Mean and standard deviation of T ratios 



s where the mean and standard deviation of the Q ratios are also known as the wild-type mean and the wild-type standard 
deviation, respectively. The mean and standard deviation of the A, C, and T means and standard deviations are also 
known collectively as the "other means and standard deviations. 

Suppose that the preceding calculations produced the foDowing data: 

m A ratios -> mean - 0.16 std. dev. - 0.003 

C ratios -> mean - 0.03 std. dev. - 0.002 

Q ratios -> mean - 0.71 std. dev. - 0.050 

16 

T ratios •> mean - 0.1 1 std. dev. - 0.004 



In one embodiment, the steps up to and including step 612 are performed In a preprocessing stage for the multiple 
wild-type experiments. The results of the preprocessing stage are stored in a f Be so that the reference calculations do 
so not have to be repeatedly calculated, improving performance 

At step 614 the highest base intensity associated with the sample sequence is checked to see if ft has sufficient 
intensity to call the unknown base. The intensity is checked by determining H the highest Intensity unknown base is 
greater than the background difference cutoff. If the intensity Is not greater than the background difference cutoff, the 
sample sequence fails to have sufficient intensity as shown at step 616. Otherwise, at step 618 the sample sequence 
25 passes by having sufficient intensity. 

At step 620 calculations are performed on the four background subtracted intensities of the sample sequence. The 
ratios of the background subtracted intensity of each base to the sum of the background subtracted intensities of the 
possible bases (ail four) are calculated, giving four ratios, one for each base. For consistency, the ratio associated with 
the reference wild-type base is called the wild-type ratio, with there being three •other" ratios. Thus, the following ratios 
so are calculated: 

A ratio -A/(A + C + G + T) 
Cratio = C/(A + C + G + T) 

55 

G ratio - G/(A+C + G + T) 



Tratio «T/(A + C + Q + T) 

40 where ratio G is the wild-type ratio and ratios A. C. and T are the "other ratios. 

Suppose the background subtracted intensities associated with the sample are as follows: 

A->310 

C->50 

G->26 
45 T->100 

Then, the corresponding ratios would be as follows: 



A ratio -310/(310 + 50 + 26+ 100) -0.64 
SO C ratio -50/ (310 + 50 + 26 + 100) -0.10 

G ratio -26/ (310 + 50 + 26 + 100) -0.05 



Tratio -100/(310 + 50 + 26+ 100) -051 

65 

At step 622 ff either the reference experiments or the sample sequence foiled to have sufficient intensity, the 

unknown jxiseis assigned the code N (insufficient intensity) as shown at step 624. 

At step 626 the wild-type and "other ratios associated with the sample sequence are compared to statistical expres- 
sions. The statistical expressions include tour predetermined standard deviation cutoffs, one associated with each base. 
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Thus, there is a standard deviation cutoff for each of the bases A, C. G, and T. The localized standard deviation cutoffs 
allow the unknown base to be called with higher precision because each standard deviation cutoff can be set to a different 
value. Suppose the standard deviation cutoffs are set as follows: 

A standard deviation cutoff -> 4 



C standard deviation cutoff •> 2 

Q standard deviation cutoff •> 8 

T standard deviation cutoff •> 4 

The wild-type base ratio associated with the sample is compared to a corresponding statistical expression: 

WT ratio * WT mean - (WTstd. dev. * WTbase std. dev. cutoff) 

where the WT base std. dev. cutoff is the standard deviation cutoff for the wild-type base. As the wild-type base is Q, 
the above comparison solves to the following: 

0.05 a 0.71 -(0.050* 8) 



0.05 2 0.31 



which is not a true expression (0.05 is not greater than 0.31). 

Each of the "other" ratios associated with the sample is compared to a corresponding statistical expression. 

Other ratio > Other mean + (Other std. dev. * Other base std. dev. cutoff) 

where the Other base std. dev. cutoff is the standard deviation cutoff tor the particular "other base. Thus, the above 
comparison solves to the following three expressions: 

A-> 0.64 > 0.16 + (0.003* 4) 



0.64 > 0.17 



C-> 0.10 > 0.03 + (0.002* 2) 



0.10 > 0.03 



T->0.21 >0.11+(0.004*4) 



0.21 > 0.13 



base is assigned the code of the reference wild-type base as shown at step 630. ^r^M^^^i 
At step 632 if one or more of the "other" ratios of the sample sequence were greater than their respective statistical 
expressions, the unknown base is assigned an ambiguity code that indicates the unknown base may be any one of the 
complements of these bases, including the reference wild-type. In this example^ the "other" Ratios tor A. C and Tj -we 
all greater than their corresponding statistical expression. Thus, the unknown base would be called t^o^emerrte 
of these bases, represented by the subset X Q, and A. Thus, the unknown base would be assigned the code D (meaning 
"not C") 

If none of the rafios are greater than their respective statistical expressions, the unknown base is assigned the code 
XOnsuflicient discrimination) as shown at step 638. r ^ armKtm 

The statistical method provides accurate base caning because it utilizes statistical data from multiple reference 
experiments to call the unknown base. The statistical method has also been used to implement confidence estimates 
and calling of mixed sequences. 
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The present Invention provides poonng processing which Is a method of processing reference and sample nucleic 
acid sequences together to reduce variations across individual experiments. Inthe representative embodiment discussed 
herein, the reference and sample nucleic add sequences are labeled with different fluorescent markers emitting light at 
different wavelengths. However, the nucleic acids may be labeled with other types of markers Including distinguishable 
radioactive markers. 

After the reference and sample nucleic add sequences are labeled with different color fluorescent markers, the 
labeled reference and sample nucleic acid sequences are then combined and processed together. An apparatus tor 
detecting targets labeled with different markers is provided in US. Application No. 08/1 95,889 and Is hereby incorporated 
by reference for an purposes. 

Fig. 13 illustrates the pooling processing of a reference and sample nucleic acid sequenca At step 702 a reference 
nucleic acid sequence is marked with a fluorescent dye, 6uch as fluorescein. At step 704 a sample nucleic add sequence 
Is marked with a dye that, upon excitation, emits light of a different wavelength than that of the fluorescent dye of the 
reference sequence. For example, the sample nucleic acid sequence may be marked with rhodamina Alternatively, the 
sample nudeic add sequence may be marked by attaching biotin to the sample sequence which will subsequently bind 
to streptavidin labeled with phycoerythrin. Of course, either sequence may be marked with these or other dyes or other 
Wnds of markers (e.g., radioactive) as long as the other sequence is marked with a marker that is distinguishable. 

At step 706 the labeled reference sequence and the labeled sample sequence are combined. After this step, process- 
ing continues In the same manner as for only one labeled sequenca At step 708 the sequences are fragmented. The 
fragmented nucleic add sequences are then hybridized on a chip containing probes as shown at step 710. 

At step 712 a scanner generates image files that indicate the locations where the labeled nucleic acids bound to 
the ctyp. There is typically some overlap between the two slgnala This Is corrected tor prior to further analysis, La , after 
correction, the data files correspond to "reference" and "sample." In general, the scanner generates an Image file by 
focusing exdtation light on the hybridized chip and detecting the fluorescent light that Is emitted. The marker emitting 
the fluorescent light can be identified by the wavelength of the light. For example, the fluorescence peak of fluorescein 
is about 530 nm while that of a typical rhodamine dye is about580nm 

The scanner creates an image file for the data associated with each fluorescent marker. Indicating the locations 
where the correspondingly labeled nudeic acid bound to the chip. Based upon an analysis of the fluorescence Intensities 
and locations, it becomes possible to extract Information such as the monomer sequence of DIMA or RNA. 

Fooling processing reduces variations across Individual experiments because much of the test environment Is com- 
mon. Although pooling processing has been described as being used to improve the combined processing of reference 
and sample nudeic acid sequences, the process may also be used for two reference sequences, two sample sequences, 
or multiple sequences by utilizing multiple distinguishable markers. 

Pooling processing may also be utilized with methods of the present invention of identifying mutations In a sample 
nucleic add sequence. These methods are highly accurate In Identifying single mutations, locating multiple mutations 
and removing false positives for mutations, where a false positive Is a base that has erroneously been Identified as a 
mutation. These methods utilize hybridization data from more than one base position to identify the likely position of 
mutations. The interrogation position on the probes is utilized to more accurately identify likely mutations which makes 
more efficient use of base calling methods. These methods may be advantageously combined with the base calling 
methods described herein to efficiently and accurately sequence a sample nucleic add sequenca 

As discussed earlier In reference to Fig. 8, the fluorescent intensities of cells near an interrogation position having 
a mutation are relatively dark which creates "dark regions" around the mutation. These lower fluorescent Intensities 
result because the cells at interrogation positions near a mutation do not contain probes that are perfectly complementary 
to the sample sequence. Thus, the hybridization of these probes with the sample sequence Is lower. The characteristics 
of these "dark regions" may be utilized to Identify mutations and false positives. 

For example, a sample sequence and a reference sequence were labeled with different fluorescent markers, In this 
case fluorescein and biotirvfehycoerytrwin. The sample and reference sequences are known and the sample sequence 
is identical to the reference sequence except for mutations at certain known positions. The sample and reference 
sequences were then processing together using the pooling processing described above and the sequern^ were rtybrid^ 
ized to a chip Induding wild-type probes that are perfectly complementary to the reference sequenca The chip Induded 
20-mer probes with the interrogation position of each probe being at the 12 th base position In the probe. 

Fig. 1 4A shows a graph of the scaled fluorescent intensities (photon counts) of the wild-type probes hybridizing wfth 
the sample and reference sequences. Along the bottom of the graph are numbers which represent wild-type cell positions 
on the chip. The photon counts of the probes In the wild-type cells are plotted on a logarithmic scale of 1 0 n . As shown, 
the photon counts range from 1 (representing a de minimus value) and 100,000. The photon counts tor the probes in 
the wild-type cell numbered "45" is around 1 0.000. 

At various wild-type cells, the photon count tor the probes In the cells drops to 1 or lower. For example, the photon 
counts for wild-type cells numbered 1 1. 24, 39, etc. are 1. The low photon counts are due to the fact that there are no 
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probes in these cells. The cells are left "Hank" in order to minimize diffraction edges and thus, the location of these Wank 
cells is known. Consequently, the intermittent wild-type cells that have a photon count of 1 do not represent erroneous 

^As shown in Fig. 14A. the scaled photon counts for the wild-type probes hybridizing with the sample and reference 
sequences are almost the same except for two "bubbles." A bubble 730 has a top curve defined by the photon counts 
ot the wild-type probes that hybridized with the reference sequence and a bottom curve defined by the photon counts 
of the wild-type probes that hybridized with the sample sequence. Following bubble 730, there is a section 732 where 
the photon counts tor the wild-type probes hybridizing with the sample and reference sequences are almost tnesame. 
After section 732 is another bubble 734 which again has a top curve defined by the hybridization of the reference 
sequence and the bottom curve defined by the hybridization of the sample sequence. Another partial bubble is shown 

*° ^EwhbutWein Fl^l 4A corresponds to a dark region surrounding a single mutation Because the wild-type probes 
at and surrounding a mutant position in the sample sequence contain a single base mismatch with the sample sequence, 
the hybridization is relatively tower which results In lower photon counts. Much information about the sample sequence 
may be acquired by a detailed analysis of these bubble regions. 

The width of the bubble indicates whether there is a false positive, asingle mutation or a multiple mutation, tf there 
is a single mutation, the width of the bubble should be approximately equal to the probe length. For example. Fig. 14A 
was produced utilizing 20-mer probes. Accordingly, bubbles 730 and 734 are approximately 20 wild-type cells wide 
indicating that the both these bubbles were produced by single mutations. The width of the dark region resulting from a 
single mutation is believed to be approximately equal to the probe length because each of the probes in this region have 
a single base mismatch with the sample sequence. 

If the width of the bubble Is substantially less than the probe length, the bubble may represent a false positive. For 
example, assume that at wild-type cell number 45 in Fig. 14A. the hybridization fto*^^^***^ 
sequence was very tow (e.g.. around 1000 photon counts). A base calling algorithm that calls the bases according to 
the intensities among the cells at lhat position may indicate that there is a mutation at the postoa i^wever. thetow 
photon counts may be due to dust on the chip and not due to lower hybridization. Since the width of* « b"bbte vrojld 
be 1 . which is substantially lower than the probe width of 20. the lower photon count at wild-type cell 45 would not be 
due to a mutation (i.e.. there 're no dark region surrounding that position). 

If the width of the bubble is substantially more than the probe length, the bubble may represent multiple mutations. 
In other words, the bubble may be produced by more than one overlapping dark region. The analysis of such a bubble 
will be discussed in more detail in reference to FIg.14C. _ .. . Ha/ , rao 

Returning to Fig. 14A. each of bubbles 730 and 734 are approximately 20 bases wide Indicating with a high degree 
of certainty that each of the bubbles represent a single rnutatioa Furthermore, the butt^^^ 

the probable location of the mutations within the bubbles. As mentioned earlier, the 20-mer probes on the chto had an 
interrogation position at the 12 th base position in the probe. Thus, the base at the 12* base position is *e baseth* 
variesamong the related WT-. A-. C-. G- and T-cells. Accordingly, the mutation should be located at the 1 2 position m 
ths bubble 

The artual mutation in bubble 730 occurs at the 12th position (from the left). AddHonaOy, the actual mutation In 
bubble 734 occurs at the 12th position (from the left). Thus, as the graph shows, there are 11 bases to the left of each 
mutation and 8 bases to the right of each mutation. By utilizing the location of the imerrogationpositiomMWn the probes, 
the present invention can help to identify the probable location of a mutation within a dark region or bubble. 

Additionally because this method identifies specific locations that may have a mutation, more efficient base calling 
may be achieved. For example, an analysis of bubble 730 indicates that there Is likely to be a single mutation around 
wild-type cell 15. Typically, most errors in base calling occur in the dark regions surrounding a mutation. Many false 
positives in this dark zone can now be eliminated because they are incompatible with the bubble sizer (which i ndfeates 
single mutation, for example). Ateo. by identifying clearly a "mismatch zone." we can now apply algorithms that factor In 
the effect of a mismatch or multiple mismatches. i»,onh 

Additionally, the shape of the bubble may indicate what mutation has occurred. Fig. 14B shows a hypotheticalgrapn 
of the fluorescent intensities vs. cell locations for wild-type probes hybridizing with two sample sequences andone 
reference sequence. A C-A mismatch win be ™re destabilize to probeh^^ 

the more destabilizing C-A mismatch results in a larger volumebubble. The shape of the bubble rnay be utilized to Identify 
the particular mutation by pattern matching bubbles stored in a library. 

Ra 14C shows a graph of the fluorescent intensities (photon counts) of the wild-type probes hybrriizing wrth the 
sample and reference sequences. A single bubble 750 is flanked on either side by regions 752 and 754wWch ido not 
contain a mutation The graph was produced from a chip containing 20-mer probes with an interrogation position at base 

" ^shownbubble 750 is 27 bases wide indicating that the bubble was produced from the dark regions surrounding 
more than one mutation as 27 is greater than 20 or the length of the probes. In addition to providing information flirt 
there are multiple mutations, analysis of the bubble Indicates the probable position of two of flie mutations. Because the 



21 



EPO 717113 A2 



interrogation position is at base 12 in the 20-mer probes, one of the mutations should be around 12 bases from the left 
end of the bubble while another mutations should be around 8 bases from the right end of the bubble. And in fact, there 
Is a mutation of T to C at wiW-type cell 62 which is 1 2 bases from the left of the bubble. Additionally, there is a mutation 
of A to G at wild-type cell 69 which is 8 bases from the right of the bubble. 

The third and last mutation within bubble 750 may be identified by performing base calling methods within the bubble. 
Alternatively, the mutation may be Identified by pattern matching bubbles from a library that indicate not only the number 
of mutations but also the specific location and type of mutation. 

Fig. 15 illustrates the high level flow of one embodiment of the present invention that uses the hybridization data 
from more than one base position to Identify mutations in a sample nucleic acid sequence. After probe Intensities from 
the hybridization of wild-type probes with a sample and reference sequence are measured, the system Identifies a bubble 
region at step 780. Bubble regions are identified as regions where the hybridization of the wild-type probes to the sample 
and reference sequence differ significantly. Additionally, the reference sequence should hybridize more strongly with the 
wild-type probes since the wild-type probes will be perfectly complementary to the reference sequence 

At step 782, the system compares the base width of the bubble to the probe ler^. If the bubble width Is substantiany 
less than the probe length, the bubble does not represent a mutation at step 784. The determination of how much less 
the bubble width may vary according to experiment conditions. 

At step 786, the system compares the base width of the bubble to the probe length to determine H they are approx- 
imately equal. If the bubble width is approximately equal to the probe length, the bubble represents a single base mutation 
at step 788. Again, the determination of how dose the bubble width should be to the probe length may vary according 
to experiment conditions. 

If the bubble width is substantially more than the probe length, the bubble represents multiple mutations at step 790. 
The system performs base calling at likely locations of mutations at step 792. The likely locations of mutations are 
determined by both the width of the bubble and the location of the Interrogation position on the probes. Additionally, the 
system may analyze the pattern of the bubble to determine the specific mutations and their positions by analyzing the 
pattern of the bubble. The base calling method with the present invention may be the Intensity ratio method, reference 
method, statistical method, or any other method. 

At step 794, the system produces confidences that the mutations are identified correctly. Each confidence is deter- 
mined by how closely the experimental data matched the data expected for the mutation that was called. For example, 
if the bubble width was exactly the same as the probe length and the base calling method identified a mutation at the 
interrogation position in the probes, there is a very high likelihood or probability that the mutation was Identified correctly 
The confidence may also be produced according to how closely the bubble pattern matched the pattern tor that mutation 
or mutations In the library of patterns. 

Although In a preferred embodiment, this method of identifying mutations in a sample nucleic add sequence Is 
utilized in conjunction with poofing processing in order to reduce variations, the method may be utilized without pooling 
processing. For example, the method may be utilized effectively where the variations between separate experiments is 
minimized or the data is adjusted accordingly. Therefore, this method is not limited to the embodiment discussed above 

The present invention provides methods of accurately identifying single mutations, locating multiple mutations and 
removing false positives for mutations. These methods are advantageously performed with pooling processing and utilize 
hybridization data from more than one base position to identify the likely position of mutations. The Interrogation position 
on the probes Is also utilized to more accurately identify the likely position of mutations which makes more efficient use 
of base calling methods. 

VI. Comparative Analysis fYlgwSftfM) 

The present invention provides a method of comparative analysis and visualization of multiple experiments. The 
method allows the intensity ratio, reference, and statistical methods to be run on multiple dataf iles simultaneously. Thfa 
permits different experimental conditions, sample preparations, and analysis parameters to be compared In terms of 
their effects on sequence calling. The method also provides verification and editing functions, which are essential to 
reading sequences, as well as navigation and analysis tools. 

Fig. 16 illustrates the main screen and the associated pull down menus for comparative analysis and visualization 
of multiple experiments (SEQ ID NO:8 and SEQ ID NO*). The windows shown are from an appropriately programmed 
Sun Workstation. However, the comparative analysis software may also be implemented on or ported to a personal 
computer, including IBM PCs and compatibles, or other workstation environments. A window 802 is shown having put! 
down menus for the following functions: File 804, Edit 80S, View 808, Highlight 810, and Help 812. 

The main section of the window Is divided Into a reference sequence area 814 and a sample sequence area 816. 
The reference sequence area is where known sequences are displayed and Is divided into a reference name subarea 
-&1 8and referencQbasssubarea 820. Thar eferencMan^subarea is shown with the filenames that contain the reference 
sequences. The chip wild-type is identified by the filename with the extension -.wtr where the # indicates a unit on the 
chip. The reference base subarea contains the bases of the reference sequences. A capital C 822 is displayed to the 
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right of the reference sequence that is the chip wild-type tor the current analysis. Although the chip wild-type sequence 
has associated fluorescence intensities, the other reference sequences shown below the chip wild-type may be known 
sequences that have not been tiled on the chip. These may or may not have associated fluorescence intensities. The 
reference sequences other than the chip wild-type are used tor sequence comparisons and may be In the form of simple 
ASCII text files. ^ , _ 

Sample sequence area 816 is where sample or unknown experimental sequences are displayed tor comparison 
with the reference sequences. The sample sequence area is divided into a sample name subarea 824 and sample base 
subarea 826 The sample name subarea is shown with filenames that contain the sample sequences. The filename 
extensions Indicate the method used to call the samplesequence where ".cq#" denotes the intensity ratio method. \iq#" 
denotes the reference method, and ".sqr denotes the statistical method (# indicates the unit on the c*to). The sample 
base subarea contains the bases of the sarrple sequences. The bases of the sample sequences are identified by the 
codes previously set forth which, for the most part, conform to the IUPAC standard. 

Window 802 also contains a message panel 828. When the user selects a base with an input device in the reference 
or sample base subarea. the base becomes highlighted and the pathname of the file containing the base Is displayed 
in the message panel. The base's position in the nucleic acid sequence is also displayed in the message panel. 

In pull down menu File 804, the user is able to loadf iles of experimental sequences that have been Wed and scanned 
on a chip. There is a chip wild-type associated with each experimental sequence. The chip wild-type associated wrth 
the first experimental sequence loaded is read and shown as the chip wild-type in reference sequence area 814. The 
user is also able to load files of known nucleic acid sequences as reference sequences tor companson purpose* ^As 
before, these known reference sequences may or may not have associated probe intensity data Addfconally, In this 
menu the user Is able to seve sequences that are selected on the screen into a project file that can be loaded in at a 
later time. The project file also contains any linkage of the sequences, where sequences are linked tor conipanson 
purposes. Sequences to be saved, both reference and sample, are chosen by selecting the sequence filename with an 
input device in the reference or sample name subareaa 

In pull down menu Edit 806, the user is able to link together sequences In the reference and sample sequence 
areas. After the user has selected one reference and one or more sample sequences, the sample sequences can be 
linked to the reference sequence by selecting an entry in the pull down menu. Once the sequences are linked a flnk 
number 830 is displayed next to each of sequences of related interest Each group of linked sequence is associated 
with a unique link number, so the user can easily identify which sequences are linked together bnk.ng sequences 
permits the user to more easily compare the linked sequences. The user is also able to remove and display links from 

thB mpufl down menu View 808, the user is able to display Intensity graphs for selected bases. Once a base is selected 
in the reference or sample base subareas. the user may request an intensity graph showing the hybridized probe inten- 
sities of the selected base and a delineated neighborhood of bases near the selected basa Intensity graphs may be 
displayed for one or multiple selected bases. The user is also able to prepare oon ^^ l a ^; e ^j^* te k ^ 

Fig iTittustratesanimensitygraphw^ 
ThefnenamecomaMngtoesequencedatate^ 

probes associated with a basa Each grouping of tour vertical bars on the graph, which are labeled as a . c . g , and 
T on line 906. shows the background subtracted intensities of probes having the indicated substitution base. In one 
embodiment, the called basee are shown In red. The wild-type base is shown at line 908. the called base is shown at 
line 910, end the base position is shown at line 912. In Fig. 1 7. the base selected is at position 120. as show" by arrow 
914 The wild-type base at this position is T; however, the called base is M which means the base is either Aor C (amino). 
The'user Is able to use intensity graphs to visually compare the intensities of each of the pcsslbtecalla 

5 Tl8 illustrates multiple intensity graph windows for selected bases (SEQ ID N032. SEQ ID NO^SEQID 
N034 and SEQ ID NO:35). There are three intensity graph windows 1002. 1004. and 1006 as shown. Each window 
may be associated with adttferent experiment, wherethe sequence analyzed in the experiment may be either a reference 
(» it has associated probe intensity data as in the chip wild-type) or a sample sequenca The windows ara aligned and 
a rectangular box 1008 shows the selected bases' position in each of the sequences (position 162 m Fig. 18). The 
rectangular box aids the user in identifying the selected basea 

Referring again to Fig. 1 6, in pull down menu Highlight 8 1 0. the user is able to corrpare the sequences of fences 
and samples. At least four comparisons are available to the user, including the following: sample sequences to the chip 
wild-type sequence, sample sequences to any reference sequence, sample sequences toa^ l^ reteence 
sequences, and reference sequences to the chip wild-type sequenca For exarnpte. after me user has bnked a reference 
and sample sequence, the user can compare the bases in the linked sequences. Bases in the sample sequence that 
are different from the reference sequence win then be indicated on the display device to the user (e.g., base is shown 
in a different color) . In another example, th e user is able to perform a comparison that will help identify sample sequences. 
After a sample is linked to multiple reference sequences, each base in the sample sequence that does notmatehthe 
wild-type sequence is checked to see if it matches one of the linked reference sequences. The bases that match a linked 
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reference sequence will then be indicated on the display device to the user. The user may then more easily identify the 
sample sequence as being one of the reference sequences. 

In pull down menu Help 812, the user is able to get information and instructions regarding the comparative analysis 
program, the calling methods, and the IUPAC definitions used in the program. 

.H 9 : - 9 i!'^ r !! eS the ln1ensity rati0 metnod cofectfy caning a mutation In solutions with varying concentrations (SEQ 
D N0:10. SEQ ID N0:11. SEQ ID N0:12. SEQ ID N0:13. SEQ ID N0:14. SEQ ID N0:15, SEQ ID N0:16 SEQ ID 
N0:17. and SEQ ID N0:18). A window 1 102 is shown with a chip wild-type 1104 and a mutant sequence 1 106 The 
mutant sequence differs from the chip wild-type at the position Indicated by the rectangular box 1 108. The chip wild- 
type and mutant sequences are a region of HIV Pol Gene spanning mutations occurring In AZT drug therapy. 

There are seven sample sequences that are called using the intensity ratio method. The sample sequences are 
actuary solutions of different proportions of the chip wild-type sequence and the mutant sequence. Thus, there are 
sample solutions 1 1 10. 1 1 1 2. 1 1 14. 1 1 16. 1 1 18. 1 120. and 1 122. The solutions are 15-mer tilings across the chip wild- 
type with increased percentages of the mutant sequence from 0 to 100% by weight The following shows the proportions 
of the sample solutions: 



so 
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Sample SoluUon 


Chip Wild-Type:Mutant 


1110 


100:0 


1112 


90:10 


1114 


75:25 


1116 


50:50 


1118 


25:75 


1120 


10:90 


1122 


0:100 



For example, sample solution 1114 contains 75% chip wild-type sequence and 25% mutant sequence 

Now referring to the bases called in rectangular box 1 1 08 tor the sample solutions, the intensity ratio method correctly 
calssan^esolution 1 110 as tia^ngabaseAasinthechlp^ildtypesequenca TOs to correct baausa sample solution 
1110 is 100% chip wild-type sequence. The intensity ratio method also calls sample solution 1112 as having a base A 
because the sample solution is 90% chip wild-type sequence. 

The Intensity ratio method calls the identified base in sample solutions 1 1 14 and 1 1 16 as being an R which to an 
^'! U ^ U 'L AC c ° da denotin 8 A or G (purine). This also a correct base call because the sample solution have from 
75% to 50% chifwviW type sequence and from 25% to 50% mutation sequence. Thus, the intensity ratio method correctly 
cans the base in this transition state. 

Sample solutions 1 1 18. 1 120. and 1 122 are caned by the intensity ratio method as having a mutation base Gat the 
22? i ^ tocabon - 17,16 is a correct base <*» because the sample solutions primarily consist of the mutation sequence 
{75%. 90%, and 100% respectively). Again, the Intensity ratio method correctly caned the bases. 

These experiments also show that the base calling methods of the present Invention may also be used tor solutions 
of more than one nucleic acid sequence. 

Pig. 20 illustrates the reference method correctly calling a mutant base where the intensity ratio method incorrectly 
caOed themutant base (SEQ ID NO.-36. SEQ ID NO:37. SEQ ID NO:38. and SEQ ID NO:39). There are three Intensity 
graph windows 1202. 1204. and 1206 as shown. The windows are aligned and a rectangular box 1208 outlines the 
bases ofinterest Window 1202 shows a sample sequence called using the intensity ratio method. However, the base 
in the rectangular box 1208 was incorrectly called base C. as there is actually a base A at that position The Intensity 
ratio method incorrectly called the base as C because the probe Intensity associated with base C is much higher than 
the other probe intensities. 

Window 1204 shows a reference sequence called using the Intensity ratio method. As the reference sequence Is 
known it is not necessary to know the method used to can the reference sequence. However, it is Important to have 
probe Intensities tor a reference sequence to use the reference method. The reference sequence Is called a base C at 
the position Indicated by the rectangular box 

Window 1206 shows the sample sequence called using the reference method. The reference method correctly cats 
the specified base as being base A. Thus, for some cases the reference method Is preferable to the intensity ratio method 
because it compares probe intensities of a sample sequence to probe Intensities of a reference sequence. 
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VILExflflfiteS 
Example 1 

The intensity ratio method was used in sequence analysis of various polymorphic HIV-1 clones using a protease 
were compared to results from an ABI sequencer. The results are inustrated below: 
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ABI 


Protease Chip 




Sense 


Antisense 


Sense 


Antisense 


No call 


0 


4 


9 


4 


Ambiguous 


6 


14 


17 


8 


Wrong call 


2 


3 


3 


1 


TOTAL 


8 


21 


29 
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SUMMARY 










ABI (sense) - 


99.5% 








Chip (sense) 


-98.1% 








ABI (arrtisense) - 98.6% 






Chip (antisense) -99.1% 
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Example 2 

HIV protease genotyplng was performed using the described chips and CallSeq™ irrt^ty ratio calculations. Sam- 
ples were it^ fromAIDS patients before and after ddl treatment. Results were 

P Fig. 21 illustrates the output of the ViewSeq™ program with four ^^^^^^^^f^^ 
pies (SEQ ID NO:22. SEQ ID N053. SEQ ID N024. SEQ ID NOS5. SEQ ID NO** and SEQ » N^. Notolto 
bawcnangeatposrtion207^ a mutaton 

ta8 Sr*T?eSpS?ie illustrative and not restrictive. Many variations of the inventor J 
those of stall In the arl [upon review of this disclosure. Merely by ^ «^ ^ '^2^2 
particular reference to theUluation of DNA (natural or unnaturaO. the methods can be used nthear^yasf™md,ps 
with other materials synthesized thereon, such as RNA. The scope of the invention J^^^SSS 
not with reference to the above description, but instead should be determmed wrth reference to the appended claims 
along with their full scope of equivalents. 
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SEQUENCE LISTING 



(1) GENERAL INFORMATION: 

(i) APPLICANT: 

(A) NAME: Affymax Technologies N,V. 

(B) STREET: De Ruyderkade 62 

(C) CITY: Curacao 

(E) COUNTRY: Netherlands Antilles 
(P) POSTAL CODE (ZIP) : none 

(ii) TITLE OP INVENTION: Computer -Aided Visualization and 
Analysis System for Sequence Evaluation 

" (ill)- NUMBER OF SEQUENCES: 39 

(iv) COMPUTER READABLE FORM: 

(A) MEDIUM TYPE: Floppy disk 

(B) COMPUTER: IBM PC compatible 

<C) OPERATING SYSTEM: PC - DOS /MS -DOS 

<D) SOFTWARE: Patentln Release #1,0, Version #1.25 (EPO) 



(2) INFORMATION FOR SEQ ID NO:l: 

(i) SEQUENCE CHARACTERISTICS: 

(A) LENGTH: 15 base pairs 

(B) TYPE: nucleic acid 

(C) STRAND EDNESS : single 

(D) TOPOLOGY: linear 

(ii) MOLECULE TYPE: DNA (oligonucleotide) 



(xi) SEQUENCE DESCRIPTION: SEQ ID NO:l: 
ATGTGGACAG TTGTA 
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(2) INFORMATION FOR SEQ ID NO: 2; 

(i) SEQUENCE CHARACTERISTICS: 

(A) LENGTH: 15 base pairs 

(B) TYPE: nucleic acid _ _ 

(C) STRAND EDNESS : single 

(D) TOPOLOGY: linear 

(ii) MOLECULE TYPE: DNA (oligonucleotide) 



(xi) SEQUENCE DESCRIPTION: SEQ ID NO: 2: 
ATGTGGATAG TTGTA 



(2) INFORMATION FOR SEQ ID NO: 3: 

(i) SEQUENCE CHARACTERISTICS: 

(A) LENGTH: 15 base pairs 

(B) TYPE: nucleic acid 

(C) STRANDEDNESS : single 

(D) TOPOLOGY: linear 

(ii) MOLECULE TYPE: DNA (oligonucleotide) 



(xi) SEQUENCE DESCRIPTION: SEQ ID NO: 3: 
ATGTGGAKAG TTGTA 
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(2) INFORMATION FOR SEQ ID NO: 4: 

(i) SEQUENCE CHARACTERISTICS: 

(A) LENGTH: 11 base pairs 

(B) TYPE: nucleic acid 

(C) STRANDEDNESS : single 

(D) TOPOLOGY: linear 

(ii) MOLECULE TYPE: DNA (oligonucleotide) 



(xi) SEQUENCE DESCRIPTION: SEQ ID NO: 4: 
AAAACTGAAA A 



(2) INFORMATION FOR SEQ ID NO: 5: 

(i) SEQUENCE CHARACTERISTICS: 

(A) LENGTH: 11 base pairs 

(B) TYPE: nucleic acid 

(C) STRANDEDNESS: single 

(D) TOPOLOGY: linear 

(ii) MOLECULE TYPE: DNA (oligonucleotide) 



(xi) SEQUENCE DESCRIPTION: SEQ ID NO: 5: 
AAAACCGAAA A 
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(2) INFORMATION FOR SEQ ID NO: 6: 

(i) SEQUENCE CHARACTERISTICS: 

(A) LENGTH: 17 base pairs 

(B) TYPE: nucleic acid 

(C) STRANDEDNESS: single 

(D) TOPOLOGY: linear 

(ii) MOLECULE TYPE: DNA (oligonucleotide) 



(xi) SEQUENCE DESCRIPTION: SEQ ID NO: 6: 
AAACCCAATC CACATCA 



(2) INFORMATION FOR SEQ ID NO: 7: 

(i) SEQUENCE CHARACTERISTICS: 

(A) LENGTH: 17 base pairs 

(B) TYPE: nucleic acid 

(C) STRANDEDNESS: single 

(D) TOPOLOGY: linear 

(ii) MOLECULE TYPE: DNA (oligonucleotide) 



(xi) SEQUENCE DESCRIPTION: SEQ ID NO:7: 
AAACCCAGTC CACATCA 
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(2) INFORMATION FOR SEQ ID NO: 8: 

(i) SEQUENCE CHARACTERISTICS: 

(A) LENGTH: 31 base pairs 

(B) TYPE: nucleic acid 

(C) STRANDEDNESS: single 

(D) TOPOLOGY: linear 

(ii) MOLECULE TYPE: DNA (oligonucleotide) 



(Xi) SEQUENCE DESCRIPTION: SEQ ID NO: 8: 
GGGGAAGCAG -ATTTGGGTAC CACCCAAGTA T 



(2) INFORMATION FOR SEQ ID NO: 9: 

(i) SEQUENCE CHARACTERISTICS: 

(A) LENGTH: 31 base pairs 

(B) TYPE: nucleic acid 

(C) STRANDEDNESS: single 

(D) TOPOLOGY: linear 

(ii) MOLECULE TYPE: DNA (oligonucleotide) 



(xi) SEQUENCE DESCRIPTION: SEQ ID NO: 9: 
GGGGAAGCAG ATTTGAAMAC CACCCAAGTA T 
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(2) INFORMATION FOR SEQ ID NO: 10: 

(i) SEQUENCE CHARACTERISTICS: 

(A) LENGTH: 59 base pairs 

(B) TYPE: nucleic acid 

(C) STRAND EDNESS : single 

(D) TOPOLOGY: linear 

(ii) MOLECULE TYPE: DNA (oligonucleotide) 



(Xi) SEQUENCE DESCRIPTION: SEQ ID NO: 10: 

GCATTAGTAG AGATATGTAC AGAAATGGAA AAGGAAGGGA AAATTTCAAA 
AATTGGGCC 



(2) INFORMATION FOR SEQ ID NO: 11: 

(i) SEQUENCE CHARACTERISTICS: 

(A) LENGTH: 59 base pairs 

(B) TYPE: nucleic acid 

(C) STRANDEDNESS: single 

(D) TOPOLOGY: linear 

(ii) MOLECULE TYPE: DNA (oligonucleotide) 



(xi) SEQUENCE DESCRIPTION: SEQ ID NO: 11: 

GCATTAGTAG AAATTTGTAC AGAGATGGAA AAGGAAGGGA AAATTTCAAA 
AATTGGGCC 
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(2) INFORMATION FOR SEQ ID NO: 12: 

(i) SEQUENCE CHARACTERISTICS: 

(A) LENGTH: 59 base pairs 

(B) TYPE: nucleic acid 

(C) STRANDEDNESS: single 

(D) TOPOLOGY: linear 

(ii) MOLECULE TYPE: DNA (oligonucleotide) 



(Xi) SEQUENCE DESCRIPTION: SEQ ID NO: 12: 

GCATTAGTAG AGATATGGAG AGRARDGGRA AXXXAAGGGA AAATTNNNAA 
AATTGGGCC 



(2) INFORMATION FOR SEQ ID NO: 13: 

(i) SEQUENCE CHARACTERISTICS: 

(A) LENGTH: 59 base pairs 

(B) TYPE: nucleic acid 

(C) STRANDEDNESS: single 

(D) TOPOLOGY: linear 

(ii) MOLECULE TYPE: DNA (oligonucleotide) 



(xi) SEQUENCE DESCRIPTION: SEQ ID NO: 13: 
^ATTAGTAG AGATATGKAS AGRARDGGRA AXXXAAGGGA AAAKTNNNAA 
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(2) INFORMATION FOR SEQ ID NO: 14: 

(i) SEQUENCE CHARACTERISTICS: 

(A) LENGTH: 59 base pairs 

(B) TYPE: nucleic acid 

(C) STRAND EDNESS : single 

(D) TOPOLOGY: linear 

(ii) MOLECULE TYPE: DNA (oligonucleotide) 



(xi) SEQUENCE DESCRIPTION: SEQ ID NO: 14: 

GCATTAGTAG AGATATGKAS AGRRRDGGRA AXXXAAGGGA AAADTYNNAA 
AATTGGGCC 



(2) INFORMATION FOR SEQ ID NO: 15: 

(i) SEQUENCE CHARACTERISTICS: 

(A) LENGTH: 59 base pairs 

(B) TYPE: nucleic acid 

(C) STRANDEDNESS: single 

(D) TOPOLOGY: linear 

(ii) MOLECULE TYPE: DNA (oligonucleotide) 



(xi) SEQUENCE DESCRIPTION: SEQ ID NO: 15: 

GCATTAGTAG AGATATGTAS AGRRADGGAA AXGGAAGGGA AAATTNNNNA 
AATTGGGCC 
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(2) INFORMATION FOR SEQ ID NO: 16: 

(i) SEQUENCE CHARACTERISTICS: 

(A) LENGTH: 59 base pairs 

(B) TYPE: nucleic acid 

(C) STRANDEDNESS : single 

(D) TOPOLOGY : linear 

(ii) MOLECULE TYPE: DNA (oligonucleotide) 



(Xi) SEQUENCE DESCRIPTION: SEQ ID NO: 16: 

GCATTAGTAG AGATATGTAC AGRGAGGGAA AXGGAAGGGA AAATTNNNNA 
AATTGGGCC 



(2) INFORMATION FOR SEQ ID NO: 17: 

(i) SEQUENCE CHARACTERISTICS: 

(A) LENGTH: 59 base pairs 

(B) TYPE: nucleic acid 

(C) STRANDEDNESS: single 

(D) TOPOLOGY: linear 

(ii) MOLECULE TYPE: DNA (oligonucleotide) 



(xi) SEQUENCE DESCRIPTION: SEQ ID NO: 17: 

GCATTAGTAG AGATATGTAS AGRGAGGGAA AXGGAAGGGA AAATTNNNNA 
AATTGGGCC 
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(2) INFORMATION FOR SEQ ID NO: 18: 

5 (i) SEQUENCE CHARACTERISTICS: 

(A) LENGTH: 59 base pairs 

(B) TYPE: nucleic acid 

(C) STRANDEDNESS: single 

(D) TOPOLOGY: linear 

70 

(ii) MOLECULE TYPE: DNA (oligonucleotide) 



(xi) SEQUENCE DESCRIPTION: SEQ ID NO: 18: 

GCATTAGTAG GAGGNNNGAC AGGGRKGGAA AXXMAAGGGA AAAKTNNNAA 
AATTGGGCC 



25 (2) INFORMATION FOR SEQ ID NO: 19: 

(i) SEQUENCE CHARACTERISTICS: 

(A) LENGTH: 160 base pairs 

(B) TYPE: nucleic acid 

so (C) STRANDEDNESS: single 

(D) TOPOLOGY: linear 

(ii) MOLECULE TYPE: DNA (oligonucleotide) 



(xi) SEQUENCE DESCRIPTION: SEQ ID NO: 19: 

TCGAGATAAT CTATGTCCTC GTCTACTATG TCATAATCTT CTTTACTTAA 
ACGGTCCTTT 60 

TACCTTTGGT TTTTACTATC CCCCTTAACC TCCAAAATAG TTTCATTCTG 
TCATGCTAGT 120 

CTATGGACAT CTTTAGACAC CTGTATTTCG ATATCCATGT 160 



60 
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(2) INFORMATION FOR SEQ ID NO: 20: 

(i) SEQUENCE CHARACTERISTICS: 

(A) LENGTH: 160 base pairs 

(B) TYPE: nucleic acid 

(C) STRANDEDNESS: single 

(D) TOPOLOGY : linear 

(ii) MOLECULE TYPE: DNA (oligonucleotide) 



15 (xi) SEQUENCE DESCRIPTION: SEQ ID NO: 20: 

NNGAGATANN NTATGTCCTC GTCYACTATG TNANNNNNNN NNNNNNNNAA 
ACGGTCCTNN 60 

20 NNNNNNNNNN NNNNNNNNNN CNNCNTAACC TCCAAAATAN NNNNNNTCTN 

NNNNANNNNT 120 

CTANNNGNAG NNNNAGANAR NCCNNNNNNN NNATNCATGT 160 

28 



(2) INFORMATION FOR SEQ ID NO: 21: 

(i) SEQUENCE CHARACTERISTICS: 

(A) LENGTH: 160 base pairs 

(B) TYPE: nucleic acid 

(C) STRANDEDNESS: single 

(D) TOPOLOGY: linear 

(ii) MOLECULE TYPE: DNA (oligonucleotide) 



(Xi) SEQUENCE DESCRIPTION: SEQ ID NO: 21: 

TCGAGATAAT CTATGTCCTC GTCTACTATG TCATAATNNN NNNNACTTAA 
48 ACGGTCCTTT 6 0 

TACCTTTGGT TTTTACTATC CCCCTTAACC TCCAAAATAG TTTCATTCTG 
NCATANNAGT 12 0 

80 CTATGNGNNG NNNTAGACAG NCCNNNNTCG ATATCCATGT 160 
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(2) INFORMATION FOR SEQ ID NO: 22: 

5 (i) SEQUENCE CHARACTERISTICS: 

(A) LENGTH: 160 base pairs 

(B) TYPE: nucleic acid 

(C) STRAND EDNESS : single 

(D) TOPOLOGY: linear 

10 

(ii) MOLECULE TYPE: DNA (oligonucleotide) 



T5 (xi) SEQUENCE DESCRIPTION : SEQ ID NO: 22: 

TCGAGATAAT CTATGTCCTC GTCTACTATG TCATAATCTT CTTTACTTAA 
ACGGTCCTTT 60 

20 TACCTTTGGT TTTTACTATC CNNCTTAACC TCCAAAATAG TTTCATTCTG 

TCATACTAGT 120 

CTATGGGTAG CTTTAGACCN CCGTATTTCG ATATCCATGT 160 



(2) INFORMATION FOR SEQ ID NO: 23: 

so 

(i) SEQUENCE CHARACTERISTICS: 

(A) LENGTH: 160 base pairs 

(B) TYPE: nucleic acid 

(C) STRANDEDNESS: single 
35 (D) TOPOLOGY: linear 

(ii) MOLECULE TYPE: DNA (oligonucleotide) 



(Xi) SEQUENCE DESCRIPTION: SEQ ID NO: 23: 

TCGAGATAAT CTATGTCCTC GTCTACTATG TCATAATCTT CTTTACTTAA 
ACGGTCCTTT 60 

TACCTTTGGT TTTTACTATC CCNCTTAACC TCCAAAATAG TTTCATTCTG 
TCATACTAGT 120 

CTATGGGTAG CTTTAGACCC CCGTATTTCG ATATCCATGT 160 
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(2) INFORMATION FOR SEQ ID NO: 24: 

(i) SEQUENCE CHARACTERISTICS: 

(A) LENGTH: 160 base pairs 

(B) TYPE: nucleic acid 

(C) STRANDEDNESS: single 

(D) TOPOLOGY: linear 

(ii) MOLECULE TYPE: DNA (oligonucleotide) 



16 (Xi) SEQUENCE DESCRIPTION: SEQ ID NO: 24: 

NCGGGATANT NTATGTCCTC GTCYACTATG TCANNNNNCN NNCNNNNCAA 
ACGGTCCNCC 60 

so NNNNNCNNNN NNCNNCYANG AANCYCAACC TCCAAAATAN NNNNNNTCTN 

NNNNANNNCN 120 

CTNNNNNNAG NGNNAGACAC CTGTATNNNN NTATNCAYGT 160 



(2) INFORMATION FOR SEQ ID NO: 25: 

(i) SEQUENCE CHARACTERISTICS: 

(A) LENGTH: 160 base pairs 

(B) TYPE: nucleic acid 

(C) STRANDEDNESS: single 

(D) TOPOLOGY: linear 

(ii) MOLECULE TYPE: DNA (oligonucleotide) 



(Xi) SEQUENCE DESCRIPTION: SEQ ID NO: 25: 

TCGRGATAAT CTATGTCCTC GTCTACTATG TCATAATCCN NNCNNCTCAA 
48 ACGGTCCTYC 60 

CNNNNYTGGT TNYTACTATC CCCCTTAACC TCCAAAATAG TTTCATTCTG 
NCATACNNST X20 

60 CTANNNNNAG NGTTAGACAC CTGTATTTCG ATATCCATGT 160 
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(2) INFORMATION FOR SEQ ID NO: 26: 

(i) SEQUENCE CHARACTERISTICS: 

(A) LENGTH: 160 base pairs 

(B) TYPE: nucleic acid 

(C) STRAND EDNESS : single 

(D) TOPOLOGY: linear 

(ii) MOLECULE TYPE: DNA (oligonucleotide) 



18 (xi) SEQUENCE DESCRIPTION: SEQ ID NO: 26: 

TCGAGATAAT CTATGTCCTC GTCTACTATG TCATAATCCN NCCTACTCAA 
ACGGTCCTTC 60 

20 TACCTTTGGT TTTTACTATC CMCCTTAACC TCCAAAATAG TTTCATTCTG 

TCATACTAGT 120 

CTATGAGTAG CTTTAGACAC CTGTATTTCG ATATCCATGT 160 



so (2) INFORMATION FOR SEQ ID NO: 27: 

(i) SEQUENCE CHARACTERISTICS: 

(A) LENGTH: 160 base pairs 

(B) TYPE: nucleic acid 

ss (C) STRANDEDNESS : single 

(D) TOPOLOGY: linear 

(ii) MOLECULE TYPE: DNA (oligonucleotide) 



(Xi) SEQUENCE DESCRIPTION : SEQ ID NO: 27: 

TCGAGATAAT CTATGTCCTC GTCTACTATG TCATAATCTT CTTTACYCAA 
45 ACGGTCCTXC 60 

TACCTTTGGT TTTTACTATC CCMCTTAACC TCCAAAATAG TTTCATTCTG 
TCATACTAGT 120 

w CTATGAGTAG CTTTAGACAC CTGTATTTCG ATATCCATGT 160 



55 
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(2) INFORMATION FOR SEQ ID NO: 28: 

(i) SEQUENCE CHARACTERISTICS: 

(A) LENGTH: 17 base pairs 

(B) TYPE: nucleic acid 

(C) STRANDEDNESS: single 

(D) TOPOLOGY: linear 

(ii) MOLECULE TYPE: DNA (oligonucleotide) 



(xi) SEQUENCE DESCRIPTION: SEQ ID NO: 28: 
AAACCCAATC CACATCM 



(2) INFORMATION FOR SEQ ID NO: 29 : 

(i) SEQUENCE CHARACTERISTICS: 

(A) LENGTH: 17 base pairs 

(B) TYPE: nucleic acid 

(C) STRAND EDNESS : single 

(D) TOPOLOGY: linear 

(ii) MOLECULE TYPE: DNA (oligonucleotide) 



(xi) SEQUENCE DESCRIPTION: SEQ ID NO: 29: 
MMACNCANNC CACANNM 
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(2) INFORMATION FOR SEQ ID NO: 30: 

(i) SEQUENCE CHARACTERISTICS: 

(A) LENGTH: 11 base pairs 

(B) TYPE: nucleic acid 

(C) STRANDEDNESS: single 

(D) TOPOLOGY: linear 

(ii) MOLECULE TYPE: DNA (oligonucleotide) 



(xi) SEQUENCE DESCRIPTION: SEQ ID NO: 30: 
TTGGGTACCA C 



(2) INFORMATION FOR SEQ ID NO: 31: 

(i) SEQUENCE CHARACTERISTICS: 

(A) LENGTH: 11 base pairs 

(B) TYPE: nucleic acid 

(C) STRANDEDNESS: single 

(D) TOPOLOGY: linear 

(ii) MOLECULE TYPE: DNA (oligonucleotide) 



(xi) SEQUENCE DESCRIPTION: SEQ ID NO: 31: 
TTGAAMACCA C 
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(2) INFORMATION FOR SEQ ID NO: 32: 

(i) SEQUENCE CHARACTERISTICS: 

(A) LENGTH: 11 base pairs 

(B) TYPE: nucleic acid 

(C) STRANDED NESS : single 

(D) TOPOLOGY: linear 

(ii) MOLECULE TYPE: DNA (oligonucleotide) 



(Xi) SEQUENCE DESCRIPTION: SEQ ID NO: 32: 
ACAGAAATGG A 



(2) INFORMATION FOR SEQ ID NO: 33: 

(i) SEQUENCE CHARACTERISTICS: 

(A) LENGTH: 11 base pairs 

(B) TYPE: nucleic acid 

(C) STRANDEDNESS : single 

(D) TOPOLOGY: linear 

(ii) MOLECULE TYPE: DNA (oligonucleotide) 



(Xi) SEQUENCE DESCRIPTION: SEQ ID NO: 33: 
AGAGRATDGG R 
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(2) INFORMATIOK FOR SEQ ID NO: 34: 

(i) SEQUENCE CHARACTERISTICS : 

(A) LENGTH: 11 base pairs 

(B) TYPE: nucleic acid 

(C) STRANDEDNESS: single 

(D) TOPOLOGY: linear 

(ii) MOLECULE TYPE: DNA (oligonucleotide) 



(xi) SEQUENCE DESCRIPTION: SEQ ID NO: 34: 
ASAGRRADGG A 



(2) INFORMATION FOR SEQ ID NO: 35: 

(i) SEQUENCE CHARACTERISTICS: 

(A) LENGTH: 11 base pairs 

(B) TYPE: nucleic acid 

(C) STRANDEDNESS: single 

(D) TOPOLOGY: linear 

(ii) MOLECULE TYPE: DNA (oligonucleotide) 



(Xi) SEQUENCE DESCRIPTION: SEQ ID NO: 35: 
ACAGGGRRGG A 
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(2) INFORMATION FOR SEQ ID NO: 36: 

(i) SEQUENCE CHARACTERISTICS: 

(A) LENGTH: 11 base pairs 

(B) TYPE: nucleic acid 

(C) STRANDEDNESS: single 

(D) TOPOLOGY: linear 

(ii) MOLECULE TYPE: DNA (oligonucleotide) 



(Xi) SEQUENCE DESCRIPTION: SEQ ID NO: 36: 
CTGGGGGGTA T 



(2) INFORMATION FOR SEQ ID NO: 37: 

(i) SEQUENCE CHARACTERISTICS: 

(A) LENGTH: 11 base pairs 

(B) TYPE: nucleic acid 

(C) STRANDEDNESS: single 

(D) TOPOLOGY: linear 

(ii) MOLECULE TYPE: DNA (oligonucleotide) 



(xi) SEQUENCE DESCRIPTION: SEQ ID NO: 37: 
CTGGCCSGTG T 



(2) INFORMATION FOR SEQ ID NO: 38: 

(i) SEQUENCE CHARACTERISTICS: 

(A) LENGTH: 11 base pairs 

(B) TYPE: nucleic acid 

(C) STRANDEDNESS: single 

(D) TOPOLOGY: linear 

(ii) MOLECULE TYPE: DNA (oligonucleotide) 



(Xi) SEQUENCE DESCRIPTION: SEQ ID NO: 38: 
CTGGGCGGTA T 
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(2) INFORMATION FOR SEQ ID NO: 39: 

(i) SEQUENCE CHARACTERISTICS: 

(A) LENGTH: 11 base pairs 

(B) TYPE: nucleic acid 

(C) STRANDEDNESS: single 

(D) TOPOLOGY: linear 

(ii) MOLECULE TYPE: DNA (oligonucleotide) 



(xi) SEQUENCE DESCRIPTION: SEQ ID NO: 39: 
CTGGCACGTG T 



Claims 

1. In a computer system, a method of identifying an unknown base in a sample nucleic add sequence, said method 

< ^ r Tnputt^nfl a^lurality of probe intensities, each of said probe Intensities being associated with a nudeic acid 

Pr0b6 ' said computer system comparing said plurality of probe intensities wherein each of said plurality of probe 
intensities is substantially proportional to said associated nucleic add probe hybridizing with at least one nucleic 
add sequence, said at least one nucleic acid sequence including said sample sequence; and 
calling said unknown base according to results of said comparing step. 

2. In a computer system, a method of identifying an unknown base in a sample nudeic add sequence, said method 

W ^^Sng aTluSty of probe Intensities, each of said probe intensities being assodated with a nucleic add 

Probe< computer system comparing said plurality of probe intensities wherein each of said plurality of prcbe 
intensities is substantially proportional to said associated nucleic add probe hybridizing with said sample sequence; 
and 

calling said unknown base according to results of said comparing step. 

3. The method of claim 2. wherein said comparing step includes the step of said computer system calculating a ratio 
of a higher probe Intensity to a lower probe intensity. 

4 The method of dahn 3, wherein said calling step includes the step of calling said unknown base according to said 
probe associated with said higher probe Intensity if said ratio is greater than a predetermined ratio value. 

5. The method of claim 4, wherein said predetermined ratio value is approximately 1 .2. 

6. In a computer system, a method of Identifying an unknown base in a sample nudeic acid sequence, said method 
comprising the steps of; . . . ^. 

Inputting a first set of probe intensities, each of said probe intensities in said first set being associated with 
a nudeic acid probe and substantially proportional to said associated nudeic acid probe hybridizing with a reference 

nudeic acid sequence; . . . . , ^ JLJ1 T nt nH 

inputting a second set of probe intensities, each of said probe intensities in said second set being associated 
with a nudeic acid probe and substantially proportional to said associated nudeic acid probe hybridizing with said 
sample sequence; 
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said computer system comparing at least one of said probe Intensities in said first set and at least one of 
said probe intensities In said second set; and 

calling said unknown base according to results of said comparing step. 

7. The method of claim 6, wherein said comparing step includes the steps of: 

calculating first ratios of a wild-type probe intensity to each probe Intensity of a probe hybridizing with said 
reference sequence, wherein said wild-type probe intensity is associated with a wild-type probe; and 

calculating second ratios of the highest probe intensity of a probe hybridizing with said sample sequence to 
each probe Intensity of a probe hybridizing with said sample sequence. 

& The method of claim 7, wherein said comparing step further includes the step of calculating third ratios of said ftet 
ratios to said second ratios. 

9. The method of claim 8, wherein said calling step includes the step of calling said unknown base according to said 
probe associated with a highest third ratio. 

10. The method of dalm 6, wherein said comparing step includes the step of calculating a ratio of a highest probe 
intensity in said first set to a highest intensity in said second set. 

1 1 . The method of claim 1 0, wherein said comparing step further includes the step of comparing said ratio of neighboring 
nucleic add probes. 

12. In a computer system, a method of identifying an unknown base in a sample nudeic add sequence, said method 
comprising the steps of : 

inputting statistics about a plurality of experiments, each of said experiments producing probe intendtiee 
each being associated with a nudeic acid probe and substantially proportional to said assodated nudeic acid probe 
hybridizing with a reference nuddc acid sequence; 

inputting a plurality of probe intensities, each of said plurality of probe intensities being assodated with a 
nudeic acid probe and substantially proportional to said associated nucleic acid probe hybridizing with said sample 
sequence; 

said computer system comparing at least one of said plurality of probe intensities with said statistics; and 
calling said unknown base according to results of said comparing step. 

13. The method of claim 12, further comprising the step of calculating said statistics. 

14. The method of claim 12, wherein said statistics \ndude a mean and standard deviation. 

15. A method of processing first and second nucleic acid sequences, comprising the steps of: 

providing a plurality of nudeic acid probes; 

labeling said first nudeic add sequence with a first marker; 

labeling said second nucleic add sequence with a second marker; and 

hybridizing said first and second labeled nudeic acid sequences at the same time. 

1 6. The method of daim 15, wherein said plurality of nudeic acid probes are on a chlpi 

17. The method of daim 15, further comprising the step of fragmenting said first and second nudeic acid sequences 
at the same time. 

1a The method of claim 15, further compridng the step of scanning for said first and second markers on said chia said 
first and second labeled nucleic acid sequences being on said chip. 

19. The method of daim 15, wherein said first and second markers are fluorescent markers that emit light at different 
wavelengths upon excitation. 

20. In a computer system, a method of identifying mutations in a sample nuddc add sequence, said method comprising 
the steps of: 

inputting a first set of probe intensities, each of said probe intensities in said first set being associated wfth 
a nudeic add probe and substantially proportional to said assodated nucleic acid probe hybridizing wfth a reference 
nudeic add sequence; 
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inputting a second set of probe intensities, each of said probe intensities in said second set being associated 
with a nucleic add probe and substantially proportional to said associated nucleic acid probe hybridizing with said 
sample sequence; 

said computer system comparing probe intensities in said first set and probe intensities in said second set 
to select hybridization regions where said probe intensities in said first set and said probe intensities in said second 
set differ; and 

identifying mutations according to characteristics of said selected regions. 

21. The method of claim 20, wherein said selected regions are determined by comparing probe intensities of wild-type 
probes. 

22. The method of claim 21 , wherein said wild-type probes are complementary to a portion of said reference sequence. 

23. Hie method of claim 21 . wherein said identifying step further includes the steps of: 

analyzing a size of a selected region; ^ _ 

identifying a likely position of a mutation in said selected region according to an interrogation position of said 
nucleic acid probes; and 

performing base calling at said likely position. 

24. In a computer system, a method of analyzing a plurality of sequences of bases, said plurality of sequences including 
" at least one reference sequence and at least one sample sequence, the method comprising the steps of: 

displaying said at least one reference sequence In a first area on a display device; and 
displaying said at least one sample sequence in a second area on said display device; 
whereby a user is capable of visually comparing said plurality of sequences. 

25. The method of claim 24, wherein said plurality of sequences are monomer strands of DNA or RNA. 

26. The method of claim 24, wherein said at least one reference sequence includes a crqp wild-type that has been tiled 
on a chip. 

27. The method of claim 26, wherein said chip wild-type sequence is displayed as a first sequence in said first area. 

28. The method of claim 26, further comprising the step of displaying a label in said first area to identify said chip wild- 
type sequence. 

29. The method of claim 24. wherein said at least one sample sequence has been hybridized on a chip. 

30. Trie method of claim 24, further comprising the step of indicate bases that dite^ 
sequences. 

31. The method of claim 24, further comprising the steps ot 

displaying a name associated with each of said at least one reference sequence In said first area; and 
displaying a name associated with each of said at least one sample sequence In said second area. 

32. The method of claim 24, further comprising the step of linking at least one reference sequence in said first area with 
at least one sample sequence in said second area. 

33. The method of claim 32, further comprising the step of indicating on said display device which sequences are linked. 

34. The method of claim 24, further comprising the step of indicating bases of said at least one sample sequence that 
' are not equal to a corresponding base in said at least one reference sequence. 

35. The method of claim 24, wherein said at least one reference sequence and said at least one sample sequence are 
aligned on said display device hybridization with said probes. 
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